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ANALYSIS  OF  SINGLE  AND  MULTIPLE  WAVELENGTH  RADAR  SYSTEMS 
FOR  ESTIMATING  SLANT  PATH  ATTENUATION  THROUGH  RAIN  AND  CLOUD  AT 
FREQUENCIES  GREATER  THAN  10  GHz 

I.O  SUMMARY  AND  CXINCLUSIONS 

The  aim  of  this  report  is  to  analyze  the  feasibility  of  using 
various  types  of  radars  including  dual  wavelength  configurations  for 
arriving  at  slant  path  attenuation  through  rain  and  cloud  at  frequencies 
up  to  100  GHz,  and  to  establish  optimum  performance  parameters  for  a 
given  system  which  will  best  perform  this  function.  The  results  associ- 
ated with  this  aim  may  be  used  in  developing  a radar  system  capable  of 
estimating  path  attenuation  through  rain  or  cloud  where  the  path  may  link 
a ground  station  and  an  aircraft  at  arbitrary  elevation  angles  and 
azimuths. 

Knowledge  of  the  path  attenuation  both  on  a case  by  case  as 
well  as  statistical  basis  is  extremely  useful  in  establishing  the  trans- 
mitter power  and  receiver  sensitivity  requirements  for  any  airplane- 
ground  link  system  (or  earth-satellite  system).  It  may  also  establish 
the  need  to  utilize  space  diversity  systems. 

Although  the  radar  systems  described  in  this  report  operate  at 
single  or  dual  wavelengths,  resulting  path  attenuations  may  be  deduced 
from  these  radar  data  at  variable  frequencies  above  10  GHz  where  the 
attenuation  due  to  rain  and  possibly  cloud  may  become  sizeable. 

The  major  analytical  topics  dealt  with  in  this  report  are: 

(1)  Parameters  associated  with  an  S-band  radar  to  measure  the 
rain  reflectivity  profile  along  a path  and  the  conversion 
of  such  profiles  to  path  attenuations  (Section  2.0). 

(2)  The  feasibility  of  using  a single  attenuating  wavelength 
radar  operating  at  15.7  GHz  (K^-band)  for  establishing 
rain  path  attenuation  at  arbitrary  frequencies  above  10 
GHz  (Section  3.0). 

(3)  The  feasibility  of  using  a K^-band  radar  for  establishing 
path  attenuation  through  cloud  at  arbitrary  frequencies 
(Section  4.0). 
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(4)  Multiple  wavelength  radar  systems  such  as  K -S  band 
(35-3  GHz)  or  band  (15.7-3  GHz)  for  establishing 
path  attenuation  at  arbitrary  frequencies  (Section  5.0). 

The  following  conclusions  have  been  arrived  at  from  the  analysis 
of  the  above  topics : 

(1)  Of  the  various  systems  examined,  the  most  straightforward 
is  an  S-band  system  monitoring  the  path  reflectivity  pro- 
files. The  conversion  to  path  attenuation  is  achieved 
using  empirical  relations  based  on  a measured  or  assumed 
drop  size  distribution. 

(2)  The  use  of  a single  attenuating  wavelength  radar  (e.g.,  one 
operating  at  15.7  GHz)  for  measuring  rain  path  attenuation 
leads  to  unacceptably  high  errors  and  should  be  avoided. 

(3)  A radar  system  operating  at  15.7  GHz  is  suitable  for  esti- 
mating cloud  liquid  water  content  from  which  cloud  path 
attenuation  may  be  determined  at  frequencies  above  10  GHz. 

(4)  A dual  wavelength  system  involving  a K and  S-band  (35-3 

d 

GHz)  radar  system  is  feasible  for  establishing  drop  size 
distributions  aloft  from  which  appropriate  empirical  rela- 
tions at  higher  frequencies  may  be  generated.  These 
empirical  relations  may  subsequently  be  used  in  a single 
wavelength  (S-band)  system  as  described  by  (1). 

(5)  A band  (15.7-3  GHz)  dual  wavelength  system  may  give 
rise  to  unacceptably  large  errors  and  its  use  should  be 
avoided. 

In  arriving  at  the  above  conclusions,  error  analyses  were 
performed  where  they  were  deemed  necessary  and  use  was  made  of  a typical 
reflectivity  profile  as  measured  at  Wallops  Island,  Virginia,  with  a 
high  resolution  S-band  radar. 

This  report  is  useful  as  a reference  for  future  work  in  the 
aforementioned  areas.  It  contains  many  derived  as  well  as  referenced 

> 
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tables  and  curves  and  the  formulations  are  either  systematically  derived 
or  referenced. 


2.0  S-BAND  RAQAR  PARAMETER  INVESTIGATION  FOR  ESTIMATING  ATTENUATION  DUE 
TO  RAIN 

We  examine  here  various  design  criteria  for  developing  2m  opti- 
mum S-band  radar  for  the  measurement  of  rain  parameters  from  which  slant 
path  reflectivity  profiles  are  obtained.  In  addition,  the  methods  by 
which  path  reflectivity  profiles  are  converted  to  integrated  path  attenu- 
ation are  examined. 

2. 1 The  Radar  Equation  at  S-Band 

We  start  with  the  expression  for  the  received  power  scattered 
from  rain  assumed  uniformly  distributed  throughout  the  pulse  volume 
(Probert-Jones  [1962]), 


P = 
r 


where 


1024  n in2 

c 

T 

A 
G 

A0  ,68. 


T A^  G^Ae^)(Aej^) 

r 


-0.2 

10  Q* 


(k  + k + k ) dr 
^ g P c' 


(2.1) 


n 

r 

k ,k  ,k 

g P c 


velocity  of  light  (3  x 10'’  m/sec) 
transmitted  power  (watts) 
pulsewidth  (sec) 
wavelength  (m) 
antenna  gain 

vertical  and  horizontal  beamwidths  (radians) 

transmitter  loss  factor  (£l.0).  (Corresponds  to 
losses  from  transmitter  power  measurement  point 
through  to  gain  measurement  point.) 

Receiver  Loss  factor  (£l.0).  (Corresponds  to  losses 
from  receiver  calibration  point  to  aptenna  gain 
measurement  point.) 

rain  reflectivity  (m"  ) 

range  (m) 

attenuation  coefficients  due  to  atmospheric  gas, 
precipitation,  and  clouds,  respectively  (dB/km). 


Since  we  are  scattering  from  rain  drops  which,  in  general,  have 
drop  diameters,  D,  given  by, 

D < .5  cm 

and  the  radar  wavelength  assumed  is. 


A = 10  cm 


hence. 


D < 1- 

A 20 


and  Rayleigh  scattering  is  applicable.  Hence,  the  reflectivity,  n,  is 
given  by  (Battan  [1973]), 


A 


where 


(2.3) 


(2.4) 


|Kol^  = 

2 

eind  where  m is  the  complex  refractive  index.  The  parameter  |k^|  has 
been  tabulated  for  various  frequencies  and  temperatures  by  Gunn  and  East 
[1954]. 

For  rain  at  A = 10  cm  for  temperatures  ranging  from  0“  to  20®C, 


(2.6) 


|K^r  “ 0.93 


(2.7) 


The  quantity,  Z,  is  called  the  reflectivity  factor  and  is  theoretically 
given  by. 


N(D)  dD 


(2.8) 
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where  D is  the  drop  diameter  and  N(D)dD  is  the  number  of  drops  per  unit 
volume  with  diameters  between  D and  D+dD  (i.e.,  the  drop  size  distribution). 
It  is  convenient  to  express  D in  mm  and  the  unit  volume  associated  with 

_3 

N(D)dD  in  m . Hence,  (2.5)  becomes, 

11^  |2  ^ ,_-18  f -1, 

n * — Z X 10  [m  ] 

X 

where  [Z]  = (mm)^/ro^  and  [X]  = m. 

Substituting  (2.7)  and  (2.9)  into  (2.1), 
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(2.9) 
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1.22  X 10 


T A0  Ae,  L L 
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(k  + k + k )dr 
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where  the  units  of  the  radar  parameters  are  given  in  the  mks  system  except 
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for  Z which  has  the  units  mm  /m  . Also,  r,  both  in  the  4enominator  and 
exponent,  is  given  in  km  and  k in  the  exponent  is  in  dB/Km.  The  resulting 
received  power,  P^,  is  expressed  here  in  milliwatts. 

2.2  Attenuation  Contributions  from  Atmospheric  Gases,  Precipitation, 
and  Cloud  at  S-Band 


It  is  convenient  to  examine  k , k 

g P 

as  the  respective  integrated  attenuations  as  a function  of  range 


and  k at  X » 10  cm  as  well 
c 


In  Fig.  1 are  given  plots  of  the  two  way  atmospheric  gas 
absorption  at  f ■ 3 GHz  for  various  elevation  angles.  These  curves  were 
obtained  from  Blake  [1969,  1970]  assuming  a CRPL  exponential  reference 
atmosphere  (N^  » 313)  where  the  absorption  is  due  to  ©2  and  H2O  collision 
broadening.  We  note  that  at  low  elevation  angles  the  atoospheric 
absorption  should  be  accounted  for  in  the  radar  equation  as  two  way 
attenuations  as  much  as  1 to  3 dB  may  result  at  the  larger  ranges. 

In  Fig.  2 is  the  attenuation  coefficient  versus  rain  rate  at 
f = 3 GHz  assuming  the  Laws  and  Parsons  drop  size  distribution  (Medhurst 
[1965]),  and  in  Fig.  3 are  plotted  the  corresponding  two  way  path  attenu- 
ations at  variable  elevation  angles  assuming  a uniform  rain  rate  up  to  an 
altitude  of  4 km  where  a 4/3  earth  radius  is  assumed  here  to  account  for 


* 


Two  way  path  attenuation,  (dB) 


Figure  1.  Two  way  attenuation  due  to  atmospheric  absorption  at  f = 3 GHz 
for  various  elevation  angles.  CRPL  exponential  reference 
atmosphere  assumed  (Ng  = 313).  Absorption  due  to  O2  and  H2O 
collision  broadening.  Normalized  to  7.5  gm/m3  at  zero  altitude 
(Blake  [1969,1970],  Skolnik  [1970]). 


Attenuation  coefficient,  (dB/km) 


T wo  way  path 


Rain  rate  (mm/hr) 


Figure  3.  Two  way  path  attenuation  versus  rain  rate  at  f = 3 GHz  at 

variable  elevation  angles  assuming  a uniform  rain  with  height 
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refraction.  It  is  apparent  that  at  low  elevations  (0®  to  2")  the  two  way 
path  attenuation  could  exceed  1 dB  for  rain  rates  greater  than  7 to  20 
mm/hr  along  the  entire  path.  However,  since  these  rain  rates  over  such 
long  ranges  are  unlikely  to  occur  even  for  widespread  rain  conditions , we 
neglect  the  attenuation  contribution  due  to  rain  at  f « 3 GHz. 

The  attenuation  contribution  due  to  water  cloud  at  X = 10  cm  is 

given  in  Table  2.1.  Assuming  the  worst  case  situation  of  0°  elevation, 

the  path  length  interval  in  a widespread  cloud  extending  from  3 to  5 km 

in  altitude  (assuming  a 4/3  earth  radius  model)  is  65.7  km.  This  gives 

_2 

rise  to  a 0.5  dB  two  way  path  attenuation  assuming  k/M  = 0.9  x 10  and 

3 

M = 0.5  gm/ra  . Since  this  attenuation  is  minimal  and  the  parameters 
selected  are  extreme  for  the  widespread  cloud  conditions,  we  neglect  the 
cloud  attenuation  at  S-band. 


Hence,  neglecting  k 


P 


and  k^  in  (2.10),  we  obtain  for  X = 10  cm. 


P 

r 


1.22  X 10"^^ 


T G LQ  A0,  L 
t e h t 


L 

g 


(2.11) 


where  the  attenuation  contribution  due  to  atmospheric  gas. 


is  given  by. 


L 

g 


(2.12) 


2.3  Dynamic  Range  Requirements 

7 6 3 

Since  reflectivity  factors  as  high  as  10  mm  /m  (i.e.,  70  dB2) 
have  been  measured  with  radar  (Donaldson  and  Tear  [1963],  Katz  [1976]), 
the  radar  dynamic  range  should  be  large  enough  to  measure  these  high 
reflectivity  targets  close  in  as  well  as  smaller  levels  (e.g.,  20-30  dBZ) 
at  greater  distances. 


ratio. 


Considering  (2.11)  at  ranges  r^^  and  r2,  we  arrive  at  the  dB 


P^Crj) 
^°«10  P^(r2) 


(2.13) 
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where  we  have  neglected  the  effect  of  gas  absorption  in  (2.13).  We  may 
define  the  left  hand  side  of  (2.13)  as  the  required  dynamic  range  of  the 
receiver  where  Z(rj^)  and  Z(r2)  represent  the  largest  and  smallest  reflec- 
tivities measured  at  the  shortest  range,  r^,  and  largest  range,  r2, 
respectively. 

In  Fig.  4 is  plotted  the  required  receiver  dynamic  range  as  a 

function  of  the  range  ratio,  r_/r  , where  Z(r,)  and  Z(r_)  are  assumed  to 
7 263  ^ i 1 2, 

be  10  and  10  mm  /m  , respectively.  Hence,  if  r^^  = 5 km,  the  required 

dynamic  ranges  at  50  and  200  km  are  70  and  82  dBZ,  respectively. 

The  above  results  suggests  the  use  of  a logarithmic  receiver, 
ideally  having  a dynamic  range  of  the  order  of  90  dB. 

2.4  Integration  Requirements 

In  arriving  at  the  integration  requirements  associated  with  the 
backscattered  signals  from  rain,  two  essential  features  of  rain  must  be 
recognized.  They  are: 

(1)  Rain  represents  a distributed  target  in  a pulse  volume, 
where  the  particles  backscatter  randomly  in  phase.  As 
such,  in  order  to  arrive  at  a reliable  estimate  of  the 
backscattered  reflectivity  a number  of  independent  samples 
must  be  averaged  for  each  pulse  volume,  Av  (i.e.,  Av  = 
(|l)(rA0^)(rA0^)). 

(2)  Because  the  rain  particles  are  in  motion  due  to  fall 
velocity  and  turbulence,  it  takes  a defined  time  interval 
for  the  particles  to  achieve  new  spatial  positions  such 
that  the  scattering  phases  have  been  sufficiently  altered 
for  causing  a decorrelated  or  independent  signal  to  be 
backscattered. 

With  regard  to  the  first  statement,  Marshall  and  Hitschfeld 
[1953]  examined  the  question  dealing  with  the  required  number  of  inde- 
pendent samples,  k,  to  achieve  a reliable  estimate  of  the  mean  backscat- 
tered power. 


Radar  dynamic  range,  (dB) 


Range  ratio  (r2/ri ) 


Figure  4.  Radar  dynamic  range  given  by  (2.13)  as  a function  of  range  ratio 
rz/ri;  ri  and  r2  are  closest  and  furthest  ranges,  respectively, 
and  Z(rjj/Z(r2)  " 
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Smith  [1966]  clarified  and  extended  some  details  of  this  work 
and  demonstrated  that  the  mean  power  can  be  arrived  at  by  selecting  only 
the  maximum  signal  from  a sample  size  of  k independent  samples.  He 
showed  this  method  to  be  more  accurate  for  k up  to  10;  however,  for 
larger  values  of  this  method  was  found  to  be  less  accurate  than  the 
straightforward  averaging  method. 

The  basic  results  of  Marshall  and  Hitschfeld  are  given  in  Figs. 
5 and  6.  Fig.  5 contains  the  probability  density  function,  P(Jj^),  of  the 
measured  average  signal  intensity,  J , where 


2 

k \ 


(2.14) 


and  represents  the  measured  k sample  of  the  intensity  (power).  The 
horizontal  scale  represents  the  measured  value  of  the  power,  Jj^.  We  note 
that  as  the  number  of  samples  averaged  reach  100  the  likelihood  increases 
significantly  that  the  measured  value  of  the  averaged  power  will  equal 

the  true  average.  Although  these  curves  were  arrived  at  analytically  for 

2 2 
A , Marshall  and  Hitschfeld  show  that  the  results  for  Log  A to  be  approx- 
imately the  same. 

Fig.  6 (derived  from  Fig.  5)  shows  the  dB  errors  in  estimating 
the  jiower  (or  log  power)  plotted  as  a function  of  sample  size  for  a 
fami iy  of  confidence  levels.  We  note,  for  example,  that  an  error  of 
10  Log  A^  within  ^0.5  dB  results  90%  of  the  time  for  a sample  size  of 
100.  For  a ssunple  size  of  20,  a 90%  confidence  level  gives  rise  to  an 
error  between  +1  and  -1.6  dB. 

With  regard  to  the  time  to  decorrelation,  theoretical  studies 
have  led  to  the  expression  (Battan  [1973],  Altas  [1964]), 

At  = 1.71  A X lO'^  sec 


(2.15) 
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Avg.  signal  intensity  (J|^) 

Figure  5.  Probability  density  function,  PCJj^)  of  the  measured 

average  signal,  Jj,.,  for  averages  comprised  of  different 
sample  sizes,  k. 


0 20  40  60  80  100 


Sample  size,  ( k) 

Figure  6.  Statistical  errors,  e , as  a function  of  sample  size, 
k,  for  different  confidence  levels. 
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where  At  is  the  time  required  for  the  autocorrelation  function  to  fall  to 
0.01  and  X is  the  radar  wavelength  in  cm.  Hence,  for  a JO  cm  radar,  it 
takes  about  20  milliseconds  after  which  the  backscattered  echoes  may  be 
said  to  be  decorrelated.  Assuming  a sample  size  of  100,  an  integration 
time  of  2 seconds  is  therefore  required  if  straightforward  methods  are 
used  in  achieving  a reliable  estimate. 

2.5  Resolution  Requirements 

2.5.1  beam  Shape  Errors 

In  measuring  a nonuniform  reflectivity  factor  with  a radar,  a 
given  inherent  smoothing  error  exists  as  the  received  echoes  correspond 
to  the  integrated  effects  of  the  true  reflectivity  integrated  throughout 
the  entire  attenna  pattern  including  sidelobes.  Donaldson  and  Tear 
[1963]  and  Donaldson  [1964a, b]  investigated  this  problem  and  examined 
errors  which  may  arise  in  both  the  estimation  of  the  reflectivities  of 
rain  cell  structures  as  well  as  errors  in  corresponding  cell  heights. 

As  an  illustration  of  the  magnitude  of  the  beam  shape  errors, 
we  show  a family  of  curves  in  Fig.  7 deduced  from  the  results  of  Donaldson 
[1964b]  for  a quadratic  circular  cell  model  defined  by, 

10  Log  Z = 50  - 9.57  x lO"^  (2.16) 

where  p is  the  radial  distance  in  kilometers  from  the  storm  center  where 
p is  confined  to  the  plane  orthogonal  to  the  radar  antenna  beam  axis. 

Although  the  original  calculations  were  made  for  model  antenna  patterns 
depicting  the  CPS-9  radar  having  1“  beamwidth,  the  results  here  have  been 
extrapolated  to  antennas  with  arbitrary  beamwidths  but  having  the  same 
general  pattern  as  used  by  Donaldson.  The  cell  model  given  by  (2.16)  has 
a peak  value  of  50  dBZ  and  has  been  scaled  from  the  original  70  dBZ  used. 

The  dashed  curve  of  Fig.  7 represents  the  true  reflectivity 
values  as  a function  of  the  distance  from  the  cell  center.  The  other 
curves  denote  the  measured  reflectivity  variations  for  the  range-beamwidth 


.law 


Measured  reflectivity  factor,  Z (dBZ) 
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0 1.0  2.0  3.0  4.0  6.0  6.0  7.0 

Distance  from  center,  p (km) 

Figure  7.  Measured  reflectivity  factor  versus  distance  from  rain  cell 
center  for  indicated  ranges  assuming  true  cell  reflectivity 
variation  shovm  (dashed  curve). 
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products  indicated.  We  note  that  for  a beamwidth  of  1“  and  a range  of 
50  nm,  a 3 dB  error  may  exist  6 km  from  the  cell  center.  This  error  is 
mitigated  by  the  fact  that  the  reflectivity  factors  here  are  small  (i.e., 
smaller  than  20  dBZ).  We,  however,  note  that  at  the  larger  ranges  and 
at  the  steeper  gradient  locations,  serious  beam  shape  errors  may  arise. 
These  results  emphasize  the  need  to  utilize  small  beamwidths  in  the 
design  of  antenna  systems  used  for  rain  reflectivity  measurements.  It 
also  demonstrates  the  limitation  of  making  measurements  at  extended 
ranges. 


2.5.2  Reflectivity  Gradient  Errors 


For  the  case  in  which  the  pulse  volume  contains  gradients  in 
the  reflectivity  factor  and  in  addition  a log  receiver  is  used,  the  proba- 
bility density  functions  depicted  in  Fig.  5 and  statistical  errors  denoted 
in  Fig.  6 are  altered.  Both  the  variance  and  the  mean  change  depending 
upon  the  reflectivity  gradient  as  well  as  the  beam  resolution.  In  Figs. 
8(a)  and  (b)  are  given  the  results  of  Rogers  [1971]  which  describe  the 
deviations  from  the  uniform  reflectivity  case  as  a function  of  Z^/Z^,  the 
ratio  of  the  power  received  from  the  strongest  reflecting  region  in  the 
pulse  volume  (proportional  to  the  reflectivity  factor,  Z^)  to  the  power 
received  from  the  weakest  portion  (proportional  to  Zj^). 

The  ordinate  parameters  in  Figs.  8(a)  and  (b)  are  defined  as, 

2 

<J>  = the  estimates  of  Log  A when  a reflectivity  gradient  is 
present. 

2 

J = the  estimate  of  Log  A for  a uniform  reflectivity  in 
pulse  volume.  This  is  equal  to  an  average  across  the 
gradient. 

2 j 

y!  * the  variance  of  the  Log  A when  a reflectivity  gradient 
^ is  present. 

2 2 

Oj  » the  variance  of  tlie  Log  A for  a uniform  reflectivity  in 
the  pulse  volume. 


<J> 


■■  Model  1 

— — — Model  2 


Reflectivity  ratio  (dB) 


iugure  8.  Deviations  of  (a)  intensity  levels  (i.e..  Log  a2)  and 
(b)  variances  of  intensity  levels  from  the  uniform 
reflectivity  case  as  a function  of  the  ratio  of  largest 
(Z2)  to  smallest  (Zj)  reflectivity  in  pulse  volume. 

The  solid  and  dashed  curves  represent  models  pertaining 
to  a linear  variation  of  Log  Z and  Z,  respectively, 
across  the  beam. 
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The  solid  aiid  dashed  curves  represent  models  pertaining  to  a linear  vari- 
ation in  log  Z and  Z,  respectively;  the  former  being  considered  as  a more 
realistic  model. 

We  again  note  the  desirability  of  utilizing  a narrow  beamwidth 

as  it  is  accompanied  by  a reduction  in  the  ratio  of  Z2/ZJ  within  the 

pulse  volume.  As  an  illustration,  it  may  be  observed  that  for  a pulse 

beam  resolution  of  1 km  and  gradients  of  reflectivity  of  10  and  20  dB/km, 

the  reflectivity  will  be  underestimated  by  1.0  and  3 dB,  respectively 

(ten  times  the  value  shown) . We  also  note  that  the  additional  variances 

2 

(Fig.  8(b)  for  these  gradients  are  0.1  and  0.33  (bels)  which  correspond 
to  standard  deviations  of  3.2  and  5.7  dB. 

2.5.3  Resolution  Errors  Due  to  Scanning 

We  assume  here  the  radar  receiver  integrator  does  not  register 

a reflectivity  value  until  a defined  integration  time,  At^,  has  elapsed. 

If  the  radar  antenna  is  scanning,  the  resolution  previously  defined  in 
part  by  the  beamwidth  alters  since  the  beam  is  sweeping  during  the  inte- 
gration time.  The  effective  resolution  may  be  defined  in  terms  of  an 

effective  beamwidth,  ^0^^,  related  to  the  true  beamwidth,  A0^,  and  scan 
rate  a [“/sec]  by  the  equation, 

A0  = A0  + aAt. 
ef  o 1 

The  error  arising  due  to  reflectivity  gradients  in  the  effective 
pulse  volume  may  be  obtained  from  Fig.  8 where  the  effective  pulse  volume 
is  defined  in  terms  of  the  effective  pulsewidth,  d0g£,  given  by  (2.17). 

2.6  Correction  Due  to  Log  Averaging  and  Processor  Quantization 

A receiver  having  a logarithmic  response  produces  an  output 

2 

power  proportional  to  the  logarithm  of  the  input  power  (i.e..  Log  A ). 

If  the  logarithm  of  the  output  power  is  sampled  and  averaged,  the 

resulting  quantity  is  the  average  of  the  logarithm  of  the  power  (i.e., 

2 

Log  rather  than  the  desired  logarithm  of  the  average  power  (i.e., 

2 

Log  A ).  Austin  and  Shaffner  [1970]  have  demonstrated  that 


(2.17) 
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10  Log  = 10  Log  + 2.5  + 


where  Q.I.  represents  the  quantization  interval  of  the  output  processor 
which  converts  the  analog  to  digital  data.  In  the  event  the  data  is 
averaged  in  analog  form  the  measured  level  underestimates  the  true 
logarithm  of  the  power  simply  by  2.5  dB.  If  a processor  is  used  which 
digitizes  the  data,  Q.I.  is  given  by, 

^ . _ dB  Dynamic  Range 

2'^  - 1 

where  N is  the  word  length  in  terms  of  the  number  of  bits  used  for  analog 

to  digital  conversion.  Hence,  a 6-bit  word  covering  a 90  dB  dynamic 

90 

range  results  in  a quantization  interval,  Q.I.  = 


2^  - 1 


5 1.4  dB. 


2. 7 Receiver  Power  Levels  from  Rain  Scattering 

We  shall  here  examine  typical  expected  received  power  levels 
for  various  conditions  of  rain  assuming  specific  radar  parameters  in  the 
bracketed  factor  of  C2.ll).  The  resulting  curves  generated  are  applicable 
however  to  arbitrary  par2uneters  by  multiplying  by  an  appropriate  factor 
(to  be  defined)  so  as  to  adjust  the  radar  constant. 

Taking  10  Log^^  of  both  sides  of  (2.11)  we  obtain  the  expression 

10  LogjQ  = dBZ  + 10  Logj^Q  K - 20  Log  r + 10  Log^Q  L^ 

where  K (or  10  Log^Q  K)  is  the  effective  radar  constant  given  by, 

K = 1.22  X 10"^^  [P^  T A0  Ae  L.  L 1 
t e r t r'' 

and  where. 


dBZ 

>■  10  Logjo 

10  Log,o  P^ 

■ receiver  power  in  dBm 

r 

a range  in  km 

10  Log^Q(L^) 

> attenuation  due  to  atmospheric  gas  (dB) 

(2.18) 


(2.19) 


(2.20) 


(2.21) 
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As  an  example  we  inject  the  following  radar  parameters  into  (2.21) 

P = lO*’  watts 

T =10  secs 


or 


Ae^  = Ae, 
e h 


= L 
t r 


Hence , 


1.75  X 10"^  rad  (1.0®) 
= 2.35  X lo"*  (43.7  dB) 

= 0.707  (1.5  dB) 


K = 1.03  X 10 


10  Log  K = -69.9  dB 


(2.22) 


The  gain  given  above  is  derived  from  the  beamwidth  assuming  a 50% 
aperture  efficiency,  e,  and  the  relations, 

Ae  = 1.2  4- 


A beamwidth  of  1®  at  A = 10  cm  hence  corresponds  to  an  antenna  diameter, 

D = 7 m. 

0 

In  Fig,  9 are  plotted  curves  describing  the  received  power 
given  by  (2.20)  as  a function  of  range  for  a family  of  reflectivity 
levels,  Z.  The  contribution  due  to  atmospheric  gas  being  at  most 
several  dB  is  here  neglected.  Although  these  curves  correspond  to  the 
above  example  parameters,  the  power  levels  may  be  converted  to  other 
system  parameters  through  the  relation. 


Py(dBm)  = P^Q(dBm)  + 69.9  + 10  Log^^  K 


(2.23) 


(2.24) 


(2.25) 


i 


9> 

5 

o 


Range,  r (km) 


Figure  9.  Received  power  (for  example  radar  system)  as  a function  of  range 
at  various  rain  reflectivity  levels.  To  convert  to  other  radar 
system  power  levels  use  (2.25). 
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where  is  the  received  power  in  dBm  for  the  given  parameters  and  K and 
P^  are  the  revised  radar  constant  and  received  power  levels,  respectively, 
for  the  alternate  system  parameters. 

In  Fig.  10  are  plotted  curves  describing  the  received  power  as 
a function  of  range  for  a family  of  effective  rain  rate  levels,  R.  We 
use  here  the  empirical  relation  [Battan,  1973], 

Z = 200  (2.26) 

which  approximately  corresponds  to  a Laws  and  Parsons  distribution  [1943] 
for  converting  the  curves  of  Fig.  9 to  those  in  Fig.  10.  In  (2.26)  Z is 
expressed  in  mm°/m  and  R is  in  mm/hr.  Here  again  the  power  levels  for 
alternate  system  parameters  may  be  obtained  using  (2.25).  We  note  from 
Fig.  10  that  a receiver  sensitivity  down  to  -100  dBm  would  enable  the 
detection  of  rain  rates  (>0.5  mm/hr)  out  to  200  km  in  range.  The  power 
levels  in  Figs.  9 and  10  may  be  adjusted  so  as  to  include  atmospheric 
attenuation  by  subtracting  the  appropriate  levels  given  in  Fig.  1. 

2.8  Estimation  of  Path  Attenuation  At  Various  Frequencies 
2.8.1  Pertinent  Regression  Relationships 

Given  a defined  drop  size  distribution,  it  may  be  demonstrated 

that  the  attenuation  coefficient,  k (in  dB/km) , may  be  related  to  the 

6 3 

reflectivity  factor,  Z (in  mm  /m  ) , through  the  empirical  relation  (Atlas 
and  U lb rich  [1974]), 

k = a Z^  (2.27) 


where  a and  b are  enpirical  parameters  dependent  strongly  on  frequency 
and  drop  size  distribution  and  less  strongly  on  temperature  as  it  influ- 
ences the  refractive  index.  Best  fit  parameters  of  a and  b may  be 
obtained  through  the  theoretical  expressions  for  Z and  k given  by, 


Z = 


N(D)dD 


[mmVm^] 


(2.28) 


Received,  power,  P,.q  (dBm) 


r “ 

-24- 


Range,  r (km) 

Figure  10.  Received  power  (for  example  radar  system)  as  a function  of  range 


at  various  rain  rate  levels.  To  convert  to  other  radar  system  ' 

power  levels  use  (2.25). 


where  Z and  the  drop  size  distribution,  N(D)dD,  have  previously  been 

3 

defined  and  C in  (2.29)  is  called  the  extinction  factor  [(dB/km)  x cm  J . 

6a  t 

It  is  a Mie  related  quantity  and  represents  a measure  of  the  attenuation 
due  to  scattering  and  absorption  from  a drop  of  diameter,  D,  for  the 
given  frequency  and  refractive  index.  These  quantities  have  been  tabu- 
lated by  several  investigators  for  the  Laws  and  Parsons  drop  size 
distribution  (e.g.,  Medliurst , (1964] ; Setzer  [1970]). 

The  form  (2.27)  may  be  arrived  at  directly  using  (2.28)  and 
(2.29)  or  indirectly  using  the  rain  rate  formulation, 

00 

R = y*  v(D)  N(D)dD  [mm/hr]  (2.30) 

o 


where  v(D)  is  the  drop  terminal  velocity  for  diameter,  D.  For  a defined 
drop  size  distribution,  empirical  relations  given  by. 


k = a R® 


(2.31) 


and 


b 

Z = a R ° 
o 


(2.32) 


may  also  be  formed  by  relating  (2.28)  and  (2.29)  with  (2.30).  The  in- 
direct method  of  arriving  at  (2.27)  is  to  combine  (2.31)  and  (2.32) 
eliminating  R.  The  parameters  in  (2.32)  are  frequency  independent  and 
only  depend  on  the  drop  size  distribution  whereas  the  parameters  in  (2.31) 
depend  on  frequency  as  well. 

A number  of  investigators  have  calculated  values  of  a and  B in 
(2.31)  for  a wide  range  of  frequencies.  More  recently,  these  parameters 


'i 


r 
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I have  been  calculated  and  tabulated  by  Rogers  and  Olsen  [1977]  and  by 

Olsen  et  al.  [1977]  for  frequencies  ranging  from  1 to  1000  GHz  for  the 
[ drop  size  distributions  of  Laws  and  Parsons  [1943],  Marshall  and  Palmer 

I [1948],  and  Joss,  Thams,  and  Waldvogel  [1968].  The  tabulations  of  Olsen 

et  al.  [1977]  in  the  frequency  range  10  through  100  GHz  are  given  in 

Table  2.2  for  the  temperatures  20°C,  0°C,  and  -10°C  (supercooled  drops). 

The  column  heads  have  the  following  significance: 

I 

r 

I LPj^  = Laws  and  Parsons  distribution  for  rain  rates  less  than 

i 50  mm/hr 

I LPjj  = Laws  and  Parsons  distribution  for  rain  rates  between 

! 25  and  150  ram/hr 

! MP  = Marshall-Palmer  distribution 

I J-T  = The  Joss  et  al.  [1968]  distribution  for  thunderstorm 

j activity 

1 J-D  = The  Joss  et  al.  distribution  for  drizzle. 

The  parameters  a and  b in  (2,32)  have  been  categorized  by 

1o  o 

Joss  et  al.  [1968]  for  drizzle  and  thunderstorm  and  are: 


140 

r1.5 

Drizzle 

500 

Thunderstorm 

Laws  and  Parsons  and  Marshall-Palmer  distributions  are  often  applied  for 
widespread  rain  conditions  and  the  corresponding  Z-R  relation  is 
approximately, 

Z “ 200  R^'^  Widespread 

If  Z.  is  considered  to  be  the  radar  reflectivity  factor  at  the 
i adjacent  range  bin  along  a radar  beam  (i^  pulse  volume)  the  esti- 
mated path  integrated  attenuation.  A,  is  given  by, 

N 

A = a Z^^  Ar 
i=l 


(2.33) 

(2.34) 

(2.35) 


(2.36) 


L 
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Table  2.2 


PARAMETER  VALUES  OF  o AND  3 IN  k = a R®.  TAKEN  FROM  OLSEN  ET  AL.  [1977]. 


T = 20 

FREQ. 

(GHz) 

a 

b 

MP 

J-T 

J-D 

MP 

J-T 

J-D 

10 

9.20x10"^ 

1. 11x10’^ 

1. 01x10'^ 

1.90x10'^ 

7.52x10’^ 

1.280 

1.230 

1.260 

1.079 

1.024 

11 

1.27x10"^ 

1.67x10’^ 

<M 

1 

O 

X 

2.42x10"^ 

9.67x10"^ 

1.259 

1.181 

1.244 

1.060 

1.047 

12 

1.68x10'^ 

2.33x10'^ 

1.81x10'^ 

3.25x10"^ 

1.22x10’^ 

1.236 

1.142 

1.223 

1.022 

1.067 

15 

3.28x10"^ 

4.59x10'^ 

3.57x10"^ 

5.89x10"^ 

2.20x10'^ 

1.173 

1.077 

1.160 

0.966 

wgm 

20 

6.83x10"^ 

8.59x10'^ 

7.51x10"^ 

0.103 

4.57x10'^ 

1.110 

1.044 

1.103 

0.934 

H9 

25 

0.113 

0.143 

0.127 

0.181 

7.76x10’^ 

1.075 

1.007 

1.064 

0.868 

1.119 

30 

0. 168 

0.228 

0.191 

0.273 

0.118 

1.043 

0.955 

1.034 

0.816 

1.113 

35 

0.235 

0.337 

0.269 

0.360 

0.167 

1.008 

0.904 

0.999 

0.782 

1.107 

40 

0.312 

0.452 

0.360 

0.439 

0.225 

0.971 

0.864 

0.966 

0.757 

1.101 

50 

0.484 

0.648 

0.572 

0.610 

0.368 

0.900 

0.816 

0.900 

0.708 

1.082 

60 

0.652 

0.775 

0.796 

0.771 

0.541 

0.845 

0.794 

0.847 

0.667 

1.055 

70 

0.802 

0.850 

1.01 

0.814 

0.735 

0.802 

0.785 

0.807 

0.663 

1.020 

80 

0.934 

0.902 

1.21 

0.810 

0.936 

0.770 

0.780 

0.772 

0.671 

90 

1.05 

0.938 

1.39 

0.876 

1.14 

0.749 

0.777 

0.744 

0.656 

0.952 

100 

1.09 

0.958 

1.53 

0.995 

1.33 

0.737 

0.774 

0.721 

0.627 

0.920 

T = 0° 

10 

1.17x10"^ 

1.14x10"^ 

1.36x10"^ 

1.69x10"^ 

1.14x10'^ 

1.178 

1.189 

1.150 

1.076 

0.968 

11 

1.50x10"^ 

1.52x10"^ 

1.73x10'^ 

2.12x10’^ 

1.41x10"^ 

1.170 

1.167 

1.143 

1.065 

0.977 

12 

1.86x10"^ 

1.96x10'^ 

2.15x10'^ 

2.62x10"^ 

1.72x10'^ 

1.162 

1.150 

1.136 

1.052 

0.985 

15 

3.21.\10"^ 

3.47x10'^ 

3.68x10'^ 

4.66x10'^ 

2.82x10"^ 

1.141 

1.119 

1.118 

1.010 

1.003 

20 

6.26x10"^ 

7.09x10'^ 

7,19x10'^ 

9.83x10'^ 

5.30x10"^ 

1.118 

1.084 

1.097 

0.946 

1.020 

25 

0.105 

0.132 

0.121 

0.173 

8.61x10“^ 

1.094 

1.029 

1.074 

0.884 

1.033 

30 

0.162 

0.226 

0.186 

0.274 

0.128 

1.060 

0.955 

1.043 

0.823 

1.044 

35 

0.232 

0.345 

0.268 

0.372 

0.180 

1.021 

0.908 

1.007 

0.783 

1.053 

40 

0.313 

0.467 

0.362 

0.451 

0.241 

0.981 

0.865 

0.972 

0.760 

1.058 
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Table  2.2 
(continued) 


T = 0°  (continued) 


FREQ. 

a 

b 

(GHz) 

MP 

1 

H 



J-D 

MP 

jgm 

J-D 

50 

0.489 

0.669 

0.579 

WEm 

0.387 

0.816 

0.905 

0.709 

1.053 

60 

0.658 

0.796 

0.801 

■Ml 

0.558 

0.850 

0.794 

0.851 

0.662 

1.035 

70 

0.801 

0.869 

1.00 

0.833 

0.740 

0.808 

0.784 

0.812 

0.661 

1.009 

80 

0.924 

0.913 

1.19 

0.809 

0.922 

0.777 

0.780 

0.781 

0.674 

90 

1.02 

0.945 

1.35 

0.857 

1.10 

0.75S 

0.776 

0.753 

0.663 

IM 

100 

1.08 

0.966 

1.48 

0.961 

1.26 

0.741 

0.774 

0.730 

0.639 

0.928 

T = -lO'C 


_ 

10 

1.32x10"^ 

1.13x10'^ 

1.59x10"^ 

1.51x10'^ 

1.42x10'^ 

1.12S 

1.171 

1.092 

1.087 

0.952 

11 

1.64x10'^ 

1.45x10’^ 

1.97x10'^ 

1.91x10’^ 

1.74x10'^ 

1.124 

1.161 

1.091 

1.077 

0.938 

12 

2.00x10"^ 

CM 

1 

O 

X 

00 

2.39x10'^ 

2.37x10'^ 

2.09x10'^ 

1.122 

1.152 

1.090 

1.066 

0.942 

15 

3.31x10'^ 

3. 14x10'" 

3.94x10'^ 

4.29x10'^ 

3.34x10“^ 

1.119 

1.135 

1.088 

1.027 

0.954 

20 

6.37x10"^ 

6.68x10'^ 

7.50x10'^ 

9.73x10'^ 

6.09x10'^ 

1.111 

1.098 

1.082 

0.951 

0.971 

25 

0. 107 

0.130 

0.126 

0.174 

9.73x10'^ 

1.091 

1.037 

1.065 

0.887 

0.988 

30 

0. 165 

0.226 

0.193 

0.276 

0.143 

1.059 

0.969 

1.037 

0.826 

1.004 

35 

0.236 

0.346 

0.275 

0.376 

0.197 

1.021 

0.910 

1.002 

0.784 

1.016 

40 

0.316 

0.468 

0.370 

0.455 

0.260 

0.981 

0.867 

0.968 

0.761 

1.025 

50 

0.487 

0.665 

0.579 

0.619 

0.404 

0.910 

0.818 

0.904 

0.714 

1.027 

60 

0.645 

0.  789 

0.787 

0.800 

0.564 

0.856 

0.796 

0.854 

0.664 

1.016 

70 

0.  777 

0.860 

0.973 

0.840 

0.728 

0.817 

0.785 

0.817 

0.659 

0.997 

80 

0.888 

0.903 

1.14 

0.808 

0.888 

0.787 

0.781 

0.789 

0.673 

0.974 

90 

0.980 

0.935 

1.28 

0.834 

1.04 

0.76S 

0.777 

0.764 

0.668 

0.953 

100 

1.05 

0.957 

1.41 

0.920 

1.18 

0.749 

0.774 

0.741 

0.646 

0.933 
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where  N is  the  nujnber  of  adjacent  pulse  volumes  in  the  supimation  and  Ar  = 

CT 

Y~  is  the  range  resolution. 

2.8.2  Validity  of  Using  Radar  for  Attenuation  Prediction 

Several  investigators  have  examined  the  validity  of  using  radar 
for  estimating  the  attenuations  at  various  frequencies  along  a slant  path 
applying  a form  similar  to  (2.36).  McCormick  [1972]  and  Strickland  [1972, 
1974]  correlated  their  3 GHz  radar  measurements  with  direct  path  measure- 
ments at  f = 15.3  GHz.  With  receiving  terminals  on  the  ground,  McCormick 
located  the  transmitter  beacon  on  an  aircraft  and  in  Strickland's  case  a 
beacon  was  located  on  the  ATS-5  satellite.  Both  investigators  reported 
good  general  correlation  between  measured  and  predicted  fade  levels, 
although  there  were  occasions  when  poor  correlations  were  noted  due  to 
possible  mixed  precipitation  phase  conditions. 

Goldhirsh  [1976a]  correlated  his  3 GHz  radar  data  with  13  and  18 
GHz  path  measurements  and  obtained  generally  good  results.  In  this  latter 
case,  13  and  18  GHz  ground  transmitters  operated  in  the  uplink  mode 
toward  the  ATS-6  satellite  and  these  reflectivities  along  the  path  were 
monitored  with  a high  resolution  radar  at  Wallops  Island,  Virginia.  At 
the  time  of  this  writing,  radar  predicted  path  attenuations  are  being 
compared  by  Goldhirsh  with  measured  attenuations  (at  f = 28.5  GHz)  using 
the  COMSTAR  satellite  beacon  2md  a ground  receiver  near  the  radar. 

Since  the  melting  layer  and  the  ice  above  it  (for  a widespread 
rain  condition)  should  have  small  influence  on  the  real  attenuation,  a 
method  was  described  by  Goldhirsh  [1976a]  for  culling  out  the  false  atten- 
uation as  predicted  by  a radar.  This  method  essentially  suggests  that 
(2.36)  be  summed  only  up  to  the  melting  layer  when  the  RHI  signature 
indicates  an  intense  layered  region  at  the  O^C  height.  When  the  RJil  shows 
a convective  rain  in  the  form  of  vertical  columns  through  the  0®C  level, 
the  summation  may  be  implemented  through  the  path. 

Using  the  form  (2.36)  with  radar  data  bases  of  rain,  several 
investigators  employed  the  use  of  radar  to  establish  fade  statistics 


i 
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(Strickland  [1974],  Zawadski  and  Rogers  [1972],  Rogers  [1972],  Inkster 
and  Rogers  [1974],  Goldhirsh  and  Robison  [1975],  Goldhirsh  [1975b], 
Goldhirsh  [1976b]  and  Hodge  [1976]).  Radar  has  the  advantage  that  it  can 
measure  the  rain  reflectivity  over  three-dimensional  space  and  using 
modeling  procedures  these  fade  statistics  may  be  converted  to  variable 
frequencies,  path  elevation  angles,  and  path  spacings  (Goldhirsh  [1975, 
1976]).  In  addition,  the  influence  of  variable  path  azimuth  and  base 
line  orientation  may  be  examined  (Goldhirsh  [1976b],  Hodge  [1976]). 

3.0  EVALUATION  OF  THE  FEASIBILITY  OF  USING  A SINGLE  ATTENUATING  WAVE- 
LENGTH RADAR  FOR  RAIN  PARAMETER  INVESTIGATION 

We  examine  here  the  difficulties  and  ambiguities  in  estimating 
rain  parameters,  using  a radar  operating  at  attenuation  wavelengths.  As 
a test  case,  a radar  operating  at  a frequency  of  15.7  GHz  is  considered. 

3. 1  Iteration  Method 

3.1.1  Derivation  of  the  Formulation 

We  derive  here  a method  for  estimating  the  rain  rate  from  the 
scatter  power  returns  of  a radar  operating  at  an  attenuating  wavelength. 
Although  the  analysis  is  performed  in  terms  of  rain  rate,  R,  it  could 
just  as  well  be  pursued  in  terms  of  the  reflectivity  factor,  Z. 

Neglecting  the  attenuation  due  to  cloud  and  atmospheric  gases, 
the  radar  equation  has  the  form  of  (2.1)  which  may  be  simply  expressed  as. 


P(r.) 


n.  -0. 

-2-^  10 

r . 
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(rj)dr 


(3.1) 


where  j denotes  the  range  bin  and  is  constant  containing  the 
multiplying  parameters  in  (2.1).  All  other  parameters  have  previously 
been  defined.  As  an  approximation,  the  Mie  parameters,  q and  k , may  be 
expressed  as. 
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b'  -1 

n = a'  R®  [m  ] 


(3.2) 


k = a 
P 


[dB/km] 


(3.3) 


where  R is  the  effective  rain  rate  in  mm/hr  and  a' , b' , d and  B are  the 
appropriately  designated  parameters. 

Substituting  (3.2)  and  (3.3)  into  (3.1)  results  in. 


P(r.)  r.^  . , -0.2  / J a[R(r)]®  dr 

[Rj]*’  10 


(3.4) 


obtain. 


Considering  the  form  (3.4)  at  the  adjacent  range  bin, 


0 ^ 


T . 

.2  f a[R(r)]^  dr 


(3.5) 


Dividing  (3.4)  by  (3.5)  and  taking  10  Logj^^  of  both  sides  of 
this  ratio,  we  obtain. 


W.  . , = U.  . , + V.  . , 
J.J  + 1 J.J  + 1 J.J  + 1 


(3.6) 


where 
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In  the  integrand  of  (3.9)  it  is  reasonable  to  assume  an  exponential  fit 
for  rain  rate  between  bins  at  r^  and  That  is, 

R(r)  = Rj  exp[-  u(r  - r^)] 


where, 


u = T-  An 
A 


R(r.) 

R(r.,i) 


A = r.  , - r. 
3 + 1 3 


Substituting  (3,10)  into  (3.9)  and  performing  the  integration. 


2 a R. 


V 

3,3  + 1 


jl  - exp[-  BuA]j 


It  is  clear  that  a knowledge  of  R^  and  R^^^  allows  for  a compu- 
tation of  Uj  as  given  by  (3.8)  and  Vj  as  given  by  (3.13). 

We  may  implement  the  iteration  technique  in  the  following  way. 

Given  a measurement  of  W.  . , and  a knowledge  of  R.  we  first  find  R.  , • 

3 .3+  A 3 3+-^ 

Once  Rj+j  Is  established,  we  repeat  the  process  by  using  the  measured 

value  of  W.  , . and  arrive  at  R-  In  this  way  a rain  rate  profile 

along  the  radar  beam  is  established.  The  iteration  technique  may  be 

pursued  using  the  following  steps  in  conjunction  with  the  relation  (3.6). 


(3.10) 
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(1)  Using  a range  bin  close  to  the  radar,  compute  the  corre- 
sponding rain  rate  from  (3.4)  where  the  attenuation  may  be  neglected 
because  of  the  short  path  length.  For  simplicity  call  this  range  bin  1. 
ilence  from  (3.4), 


‘^l  = 


\Tl  P(r^) 


1_ 

\b' 


(2)  With  and  ^ [given  by  (3.14)  and  (3.7),  respectively] 
obtained  from  the  measured  power,  assume  for  the  first  approximation 

2 = 0 3nd  compute  R2  from  ^ as  given  by  (3.8).  Call  this  rain  rate 
Ri  to  denote  the  first  approximation, 

(3)  Compute  the  first  approximation  to  ..  as  given  by  (3.13) 

2 i I i 

by  substituting  R and  R«  into  this  expression.  Call  this  first  approxi- 

1 ^ ^ 

^ ^ If* 


mation  V 


1,2* 


That  is. 


(4)  Solve  for  the  second  approximation  to  U by  using  (3.6). 

^ 9 ^ 


2 1 
U = W - V 

^1,2  ”1,2  ^1,2 

With  the  right  hand  side  known,  a second  approximation  to  R-  may  be  cal- 

2 ^ 
culated  using  (3.8);  say  R^. 

2 

(5)  Repeat  steps  (3)  and  (4)  using  the  value  of  R_  and  arrive 

3 / 

at  the  third  approximation  to  rain  rate  at  r2:  namely  R2,  This  process 
may  be  continued  N times  until  convergence;  that  is  until, 

«2  - "2''  ‘ ' 

where  e is  a predesignated  uncertainty. 

(6)  After  R2  has  been  calculated  the  rain  rate  R^  may  be 
obtained  using  steps  (2)  through  (5)  by  replacing  subscripts  1 and  2 in 


(3.14) 


(3.15) 


(3.16) 
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these  steps  by  2 and  3,  respectively.  In  this  manner,  the  entire  rain 
rate  profile  may  be  calculated. 

3.1.2  Testing  the  Method 

As  an  illustration  of  the  iteration  method,  we  select  here  a 
segment  of  rain  rate  profile  as  measured  at  Wallops  Island  with  an  S-band 
radar  (SPANDAR)  on  July  10,  1976,  at  19:52:20  GMT  at  60“  azimuth  and  20“ 
elevation.  The  rain  rates  were  computed  from  the  measured  values  of  Z 
using  the  empirical  form, 

^ = 200  (3.17) 

The  rain  rate  profile  over  an  approximate  2 km  segment  is  shown  plotted 
(solid  line)  in  Fig.  11  where  a level  is  indicated  every  0.15  km  which 
corresponds  to  the  range  resolution.  A,  of  the  radar  (1  iisec  pulsewidth) . 

Using  the  results  of  Haddock  [1956]  and  assuming  a Laws  and 
Parsons  distribution,  we  obtain  for  f = 15.7  GHz, 


a'  = 9.808  X 10"  (3.18) 

b'  = 1.438  (3.19) 

From  the  results  of  Medhurst  [1965], 

_2 

a = 4.343  X 10  (3.20) 

B = 1.122  (3.21) 


The  parameters  given  by  (3.20)  and  (3.21)  may  be  substituted  in  (3.13) 
where  A = 0.15  km. 

We  assume  the  profile  given  in  Fig.  11  (solid  line)  to  be  the 
true  rain  rates  and  compute  the  values  of  W.  . , that  would  be  measured 
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Figure  11.  Segment  of  rain  rate  profile  (solid  line)  as  measured  at  Wallops 
Island,  Va.,  on  July  10,  1976,  at  19:52:20  GMT  at  60“  azimuth 
and  20“  elevation  with  an  S-bauid  radar.  Other  points  represent 
apparent  measured  rain  rates  for  given  calibration  errors  using 
a 15.7  GHz  radar  and  iteration  method. 
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by  a 15.7  GHz  radar.  This  is  achieved  by  computing  individually  the  true 

values  of  U.  . , given  by  (3.8)  and  V.  . . given  by  (5.13).  Hence  the 
J » J*!  J « J 

computed  values  of  W.  . . are  considered  true  and  the  same  as  would  be 

J » J'*’* 

measured  with  the  15.7  GHz  radar. 


The  measured  rain  rates  at  ranges  less  than  8.85  km  were  less 
than  1 mm/hr  and  therefore  have  been  neglected  in  Fig.  H,  The  true  rain 
rate  at  the  first  iteration  range  bin  is  therefore  3.65  mro/hr  and  this  is 
the  value  a 15.7  GHz  radar  would  measure  assuming  no  uncertainties. 

Hence  we  have  a priori  knowledge  of  ^ 1^2  • We  may  now  proceed  to 

evaluate  R2  as  well  as  rain  rates  at  subsequent  range  bins. 

Selecting  an  e = 0.1  mm/hr  in  (3.16),  the  number  of  iterations 

ranged  from  4 to  6 and  the  technique  is  found  to  work  well  over  the 

entire  range  interval  giving  a rain  rate  profile  that  essentially  coin- 
cides with  the  true  values  in  Fig.  11.  However,  if  we  now  introduce 
uncertainties  into  the  calibration  constant,  in  (3.1)  of  ^1,  ^2,  and 
j^3  dU,  the  indicated  uncertainties  (between  the  respective  symbols)  in 
Fig.  11  result. 

We  note  that  for  an  uncertainty  of  ^ 1 dB,  rain  rate  errors  of 
more  than  100%  result  after  only  2 km  of  iteration.  Hence  even  for  small 
uncertainties  the  errors  become  unacceptably  large  with  increasing  range, 
■p-ic  reason  for  this  is  due  to  the  bootstrap  operation  of  the  method  where 
•’e  rain  rate  at  any  given  range  bin  is  dependent  on  the  accuracy  of  the 

rain  rate  at  the  preceding  bin.  This  results  in  a multiplication  of 

errors  that  is  dependent  on  the  number  of  rain  bins  considered. 


3.2  Errors  Arising  by  Neglecting  Attenuation 

It  is  of  interest  to  compare  the  results  in  Fig.  11  with  those 
for  the  case  in  which  attenuation  is  entirely  neglected  at  15.7  GHz. 

The  so-called  true  rain  rate  for  the  bin  may  be  obtained 
from  (3.4)  and  is. 
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The  "apparent"  rain  rate  neglecting  attenuation  is, 


fr;  P(r  ) 
R.  = 1-4 J. 


1 

kb' 


ja  \ a' 

Comparing  the  apparent  eind  true  rain  rates  we  observe, 


R.  = R.  10 
ja  J 


- ^ c.[R(r)]®  dr 


which  may  alternately  be  expressed  by. 


0.2 


n-1  r. 


R.  = R.  10 
ja  3 


a[R(r)J®  dr 


J-i  r: 


where  n is  the  number  of  range  bins  considered. 

Substituting  the  form  (3.9)  into  the  exponent  of  (3.25), 
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where  V.  . , is  given  by  (3.13).  The  apparent  rain  rate  as  given  by 

J » J k 

(3.26)  is  now  a form  ripe  for  computation. 


In  Fig.  12  are  plotted  both  the  true  rain  rate,  R^  (solid  line), 

as  well  as  the  apparent  rain  rate,  R.  (dashed  line),  given  by  (3.26). 

J a 

We  observe  that  errors  in  excess  of  100%  result  after  only  two  kilometers 
for  the  case  in  which  the  influence  of  attenuation  of  the  measured  power 
is  neglected.  It  should  be  noted  that  the  above  errors  are  in  general 
even  larger  when  calibration  or  measurement  uncertainties  are  included. 


3. 3 Hitschfeld-Bordan  Solution 

A direct  method  for  arriving  at  the  true  rain  rate  given  rain 
attenuation  has  been  described  by  Hitschfeld  and  Bordan  [1954].  Using 
their  methods,  an  equation  of  the  form, 


-39- 


dU 


^ + 3 U f (r)  + K 3 = 0 


where 


U = R. 


f(r)  . log  R.^ 


K = iii-i  In  10 


and  where  the  true  rain  rate,  , is  given  by  (3.22)  and  the  apparent 
rair 
by. 


rain  rate,  R.  , is  defined  by  (3.23).  The  solution  to  (3.27)  is  given 
J ^ 


Rj(r)  = 


R.  (r) 
j a 


K 3 


a/B 


(V)d  V 


Hence  given  the  power  received  from  rain  scattering  at  each 

range,  r.  [i.e.,  P(r.)],  R.  as  given  by  (3.23)  may  be  calculated  and  the 
J 11“ 


true  rain  rate,  R^ (r) , given  by  (3.31)  is  established. 


The  expression  (3.31)  is  exact  assuming  the  empirical  relations 
(3.2)  and  (3.3)  are  correct.  However,  the  solution  suffers  from  the  same 
difficulty  as  the  previously  described  iteration  method;  namely,  a small 
uncertainty  in,  for  example,  the  calibration  constant  results  in  unaccept- 


ably large  errors  in  Rj (r) . As  a matter  of  fact  the  two  methods  may  be 


shown  to  be  equivalent  and  the  errors  in  R^ (r)  are  the  same  as  those  using 
the  iteration  methods  for  the  same  uncertainties. 


4.0  CHARACTERISTICS  OF  CLOUD  PROPERTIES  USING  A RADAR  AT  K^-BAND 


We  examine  here  the  power  scattered  from  cloud  using  a radar  at 
15.7  GHz  and  the  ability  of  estimating  the  cloud  liquid  water  content,  M 


(3.27) 


(3.28) 

(3.29) 

(3.30) 


(3.31) 
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3 

(gm/m  ),  from  the  received  radar  echoes.  A knowledge  of  M allows  for  the 
computation  of  the  attenuation  coefficient  and  ultimately  the  integrated 
attenuation  at  other  wavelengths.  The  radar  pareimeters  at  15.7  GHz  used 
represent  nominal  values  of  a radar  under  design  at  RADC. 

4. 1 Attenuation  Contributions  from  Atmospheric  Gases 

In  using  a radar  at  15.7  GHz  at  shallow  elevation  angles,  it  is 
essential  to  establish  the  attenuation  contribution  due  to  atmospheric 
gases.  In  this  section  we  evaluate  this  contribution  using  the  theory  of 
Van  Vleck  [1947,  1949]  which  is  summarized  by  Bean  and  Dutton  [1968]. 

The  total  attenuation  coefficient  is  essentially  due  to  the  contributions 
of  oxygen  and  water  vapor.  That  is, 

kg  = kCop  + kCH^O) 


The  contribution 

due  to  0^ 

is  given  by. 

= 0-34  /_ 
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where  X is  the  wavelength  in  cm,  P is  the  pressure  in  bars,  and  T is  the 
temperature  in  "K.  Also,  and  Av^  [cm  ^]  are  line  width  factors  given 

by. 
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r. 


The  contribution  due  to  U~0  is  comprised  of  two  components;  that 

^ n ) 

arising  from  the  absorption  line  at  1.35  cm  [k|^  q]  and  that  arising 


(2) 

from  the  absorption  bands  above  the  1.35  cm  line  [k^  q].  These  contribu- 
tions are  given  by, 


k^^^(H20) 


0,0318  /293\^'^,„_/  644 \ 
( T ) T J 


/\v, 


(X'^  - 0.74074)^  + AVj 


Av, 


(X'^  + 0.74074)^  + AVj 


and 


u(2),ur,^  0.05  /293\  / P \ IZlSy-^r, 

k (H2O)  = p ( T ) \1. 01325  / \ T / 

A 

3 

where  p is  the  water  vapor  density  (gm/m  ) , Av^  is  the  line  width  factor 
at  1.35  cm,  and  Av^  is  the  effective  line  width  of  the  absorption  bands 
above  the  1.35  cm  line.  The  line  widths  are  given  by. 


Av. 


Av, 


0.087 


(r 


J—\  /iiiV 

01325 / \ T / 


0.5 


(1  + .0046  P) 


In  Fig.  13  are  given  curves  describing  the  one  way  integrated 
path  attenuation  at  15.7  QIz  due  to  atmospheric  absorption  as  a function 
of  range  and  path  elevation.  As  a basis  for  comparison,  a similar  set  of 
curves  are  given  in  Fig.  14  for  f = 35  GHz.  These  curve*  were  calculated 
using. 


r 

y [k(02)  + k^^^CH20)  + k^^^(H20)]dr 


where  the  component  contributions  in  the  integrand  of  (4,8)  are  given  by 
(4,2),  (4,5),  and  (4,6).  Temperature  and  pressure-height  profiles  were 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


in  p u) 
ih  ih 


8. 

c 

a (s 

S oc 


^ -s 


ppotno  potno 
(gp)  uo!)enu8))e  Ab/m  buq 


n 


O 0) 

"3) 
c S 

o 

'H  C 

t{  ® 

C«  .H 

3 

C 

4>  > 

^ 0) 

W •-«  . 

cd  0)  -o 
o 
B 

^ S ^ 

Cg  HS  </) 

CXt  i/) 

4>  CO 

taO 

5 §.I2 

R>  M 
M t<1 

&o*i-i  e 

4>  O ^ 
■iJ  E 
C C M 
■H  O 
•H  m 
^ 4J  • 

rt  o r'. 

* s 

3 11 

4>  <« 

£ o 

Orta 


Figure  14.  One  way  integrated  path  attenuation  due  to  atmospheric  O2  and  II2O  absorption  as 
a function  of  range  and  elevation  angle  at  f = 35  GHz.  A ground  humidity  of 
7.5  gm/ra^  is  assumed. 
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taken  from  the  U.S.  Standard  Atmosphere,  1972,  for  July  at  midlatitudes 
(Cole  et  al.  [1965])  and  are  listed  in  Table  4.1.  A 4/3  earth  model  was 
assumed  to  account  for  ray  bending  due  to  refractive  index  changes. 

Comparing  mean  monthly  values  of  precipitable  water  to  values 
of  surface  dew  point  for  stations  around  the  country,  Reitan  [1963]  found 
that  an  appropriate  vapor  density-height  profile  may  be  given  by. 


p - exp(-  0.44  h) 


[gm/m^] 


(4.9) 


where  p is  the  surface  vapor  density.  Equation  (4.9)  was  used  in  the 

° 3 

computation  of  the  curves  in  Figs.  13  and  14  where  p =7,5  gm/m  was 

o 

chosen  as  a typical  value  of  the  vapor  density.  An  HP  9825  minicomputer 
was  used  in  the  evaluation  of  the  integral  (4.8). 

It  is  interesting  to  compare  the  above  model  with  the  measured 
results  of  Altshuler  [1977]  who  used  suntracker  data  to  obtain  the  inte- 
grated path  attenuation  through  the  atmosphere  in  the  Boston  area.  In 
Figs.  15  and  16  are  given  the  total  one  way  path  attenuation  as  a function 
of  elevation  angle  for  the  above  model  used  and  Altshuler's  results  at 

15  and  35  GHz.  Two  ground  average  vapor  densities  are  used  in  the  models; 
3 3 

namely  7.5  gm/m  and  5.0  gm/m  . We  note  that  both  ground  vapor  densities 
give  results  at  15.7  GHz  that  are  well  within  the  standard  deviation  of 
the  measured  results  (Fig.  15).  At  35  GHz  (Fig.  16)  the  7,5  gm/m^  v^or 
density  gives  results  that  are  slightly  outside  the  measured  standard 
deviation  values,  whereas  the  5,0  gm/m^  values  are  well  inside. 

4.2  Attenuation  Contribution  from  Cloud 

Hie  contribution  from  a widespread  cloud  condition  especially 
at  shallow  elevation  angles  is  examined  here.  The  attenuation  due  to 
atmospheric  gas  and  cloud  is  given  by. 


A A = 
8 c 


(k  + k )dr 
g C-' 


(4.10) 
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Table  4.1 

TEMPERATURE  AND  PRESSURE  PROFILES  FOR  JULY  AT  MID  LATITUDES. 
From  U.S.  Standard  Atmosphere,  1962  (Cole  et  al.  [1965]) 


h(km) 

TC’K) 

P(atm) 

0 

294.15 

1.01325 

1 

289.65 

0,902198 

2 

285.15 

0.801594 

3 

279.16 

0.710433 

4 

273.17 

0,628063 

5 

267.17 

0.553615 

6 

261.19 

0.486632 

7 

254. 70 

0.426404 

8 

248.22 

0.372351 

9 

241.73 

0.32402 

10 

235.25 

0.28094 

One  way  attenuation  (dB) 


'1 
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Figure  15.  Measured  one  way  atmospheric  path  attenuation  by  Altshuler 
[1977]  as  a function  of  elevation  angle  compared  with  model 
atmosphere  attenuation  at  f = 15.7  GHz. 
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Figure  16.  Measured  one  way  atmospheric  path  attenuation  by  Altshuler 
[1977]  as  a function  of  elevation  angle  con^ared  with  model 
atmosphere  attenuation  at  f = 35  GHz. 


-48- 


The  attenuation  coefficient  due  to  cloud  is  given  by, 

= K [dB/km]  (4.11) 

where  K is  a proportionality  constant  [dB/km/gm/m^]  and  M is  the  cloud 

3 ^ 

liquid  water  content  [gm/m  ]•  The  proportionality  constants  have  been 
interpolated  from  the  tabulations  of  Gunn  and  East  [1954]  and  are  listed 
in  Table  4.2  for  15.7  and  35  GHz.  In  arriving  at  the  cloud  attenuation 
contribution,  we  assume  the  following:  (1)  a 4/3  earth  model,  (2)  the 
values  of  K in  Table  4.2,  (3)  the  standard  atmosphere  temperature  profile 
in  Table  4.1,  and  (4)  a cloud  layer  between  3 and  5 km.  The  cloud  contri- 
bution was  added  to  the  atmospheric  gas  contribution  evaluated  i .ection 
4.1.  The  results  are  given  in  Figs.  17  through  24  for  both  15.7  and  35  GHz 
where  the  one  way  attenuations  are  plotted  as  a function  of  range  for 
various  liquid  water  contents,  M^,  cUid  elevation  angles.  Although  liquid 
water  contents  as  high  as  1 gm/nf^  are  considered,  average  values  of  0.1 
and  0.2  gm/m’’  are  more  likely  for  the  stratus  type  clouds  (Kantor  [1965]). 

3 

We  note  that  for  nominal  liquid  water  contents  of  0.1  gm/m  , 
the  addition  contribution  of  cloud  attenuation  is  less  than  1.5  dB  (0“ 
elevation)  at  15.7  and  7 dB  (0°  elevation)  at  35  GHz.  For  increasing 
elevation  the  cloud  attenuation  reduces. 

In  Figs.  25  and  26  we  compare  the  model  one  way  overall  attenu- 
ation (cloud  and  atmospheric  gas)  with  the  cloud  attenuation  measurements 
of  Altshuler  [1977]  at  15  and  37  GHz.  We  note  that  the  modeled  results 

are  well  within  the  standard  deviations  of  the  measured  results  assuming 

3 3 

= 7.5  gm/m  and  = 0.1  gm/m  . 

4.3  Power  Scattered  from  Cloud 


We  examine  here  the  power  scattered  from  cloud  using  a radar 
operating  at  15.7  GHz.  Starting  with  the  radar  equation  (2.10),  we 
i assume  the  following  radar  parameters  (suggested  by  RADC) ; 
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Table  4.2 

INTERPOLATED  VALUES  OF  k^/M^  FROM  GUNN  AND  EAST  [1954] 


Figure  17.  One  way  path  attenuation  as  a function  of  range  at  f = 15.7  GHz,  0“  elevation, 

for  a given  atmosphere  and  cloud  model  and  variable  cloud  liquid  water  contents. 


attenuation  as  a function  of  range  at  f = 15.7  GHz,  4°  elevation, 
atmosphere  and  cloud  model  and  variable  cloud  liquid  water  contents 
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Figure  22.  One  way  path  attenuation  as  a function  of  range  at  f = 35  GHz, 
2®  elevation,  for  a given  atmosphere  and  cloud  model  and 
variable  cloud  liquid  water  contents. 
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One  way  attenuation  (dB) 


Figure  25. 


Elevation  angle  (degrees) 


Measured  one  way  path  attenuation  through  cloud  (Altshuler 
[1977])  as  a function  of  elevation  angle  coiqpared  with  cloud 
and  model  atmosphere  attenuation  at  f » IS. 7 GHz. 
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Figure  26.  Measured  one  way  path  attenuation  through  cloud  (Altshuler 

[1977])  as  a function  of  elevation  angle  compared  with  cloud 
and  model  atmosphere  attenuation  at  f » 35  GHz. 


-60- 


X = 1.91  cm 


= 2 X 10^  watts 
T = .08  usee 
G > 60  dB 

ae  = ae.  * o.is® 

e h 

= L = 1.0 
t r 


(4.12) 


Hence  (2.10)  becomes. 


-8  Z f 

P - 3.67  X 10  ^ 10  o*'  » ® 


C4.13) 


where  P is  the  received  power  in  mw,  r is  the  range  in  km,  and  Z is  the 
cloud  reflectivity  in  mm  /m  . A value  of  |k^|  = 0.92  corresponding  to 
the  temperature  of  0®C  and  frequency  of  15.7  GHz  was  used  in  (4.13)  (Gunn 
and  bast  [1954]). 

Using  the  results  of  Atlas  and  Bartnoff  [1953]  relating  Z to 

3 

[gm/m  ] for  cloud;  namely, 

Z = 4.8  X 10"^  (4.14) 


(4.13)  becomes. 
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The  received  power  as  a function  of  range  is  plotted  in  Fig.  27 
for  various  elevation  angles  and  a cloud  liquid  water  content  of  0.2  gm/m^. 

We  note  that  for  a receiver  sensitivity  of  -145  dBm  (suggested  by  RADC) 
cloud  scatter  can  be  detected  at  elevation  angles  equal  to  or  greater 
than  2®.  For  » 0.1  gm/m  , all  the  curves  shift  down  by  6 dB.  For 
this  latter  case  and  for  the  same  receiver  sensitivity  of  -145  dBm,  cloud  , 
at  path  angles  greater  than  4°  are  detectable. 

It  would  thus  appear  that  the  15.7  GHz  radar  with  the  given 
receiver  sensitivity  (-145  dBm)  represents  an  excellent  tool  for 
measuring  cloud  liquid  water  content.  At  elevation  angles  greater  than 
2®  and  4®,  the  two  way  total  attenuation  should  be  less  than  6 and  4 dB 
(assuming  the  given  cloud  model  and  = 0.2  gm/m  ) [Figs.  18  and  19]. 

The  attenuation  due  to  cloud  itself  is  observed  to  be  less  than  4 and  2 dB. 

4.4  Comparison  of  Scatter  from  Cloud  Due  to  Refractive  Index 

Variations 

We  examine  here  the  possibility  of  confusion  arising  from 
scatter  from  irregularities  of  refractive  index  (i.e.,  clear  air  echoes). 

Under  given  assumptions  of  isotoropy,  homogeneity,  and  eddy  size  dimen- 
sions, the  reflectivity  is  given  by  (Tatarski  [1964]), 

hcA  = 39  Cn^  [m'^]  (4.17) 

where 

2 -2/3 

C 3 structure  constant  [cm  ] 
n 

X » wavelength  [cm] 

The  reflectivity  due  to  cloud  scatter  at  f ■ 15.7  GHz  corre- 
sponds to  the  Rayleigh  expression. 
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Injecting  |k^|  = 0.92  and  (4.14)  into  (4.18)  results  in, 

rip,  = 1.35  X 10“^  M ^ [m“^] 

UL  C 

and  dividing  (4.19)  by  (4.17)  gives. 


]CL 

^CA 


3.47  X 10 


-11 


M 


c 2 

n 


At  f = 15.7  GHz  (X  * 1.91  cm),  (4.20)  becomes. 


-12 

3.23  X 10 


In  Fig.  28,  the  ratio  (4.21)  is  plotted  as  a function  of  the 

2 

structure  constant  for  various  values  of  M . We  note  that  even  for 

3 2 -14 

small  values  of  M >0.1  emd  0.2  gm/m  and  large  values  of  C =10 

-2/3  ^ ^ 

cm  (Hardy  and  Katz  [1969])  the  reflectivity  due  to  cloud  exceeds  that 

due  to  refractive  index  variations  by  a factor  greater  than  3 to  13, 

respectively.  Hence  little  confusion  should  arise  from  the  simultaneous 

existance  of  irregularities  of  refractive  index  imbedded  in  the  cloud 

medium. 


4.5  Extrapolation  of  Cloud  Attenuation  to  Other  Wavelengths 

Assuming  an  estimate  of  the  cloud  liquid  water  content,  is 
arrived  at  using  the  15.7  GHz  radar,  the  cloud  attenuation  may  be  obtained 
at  other  wavelengths  using, 

= K [dB/km] 

where  K may  be  extracted  from  the  plots  of  K versus  frequency  as  given 
in  Fig.  29  for  various  temperatures.  The  plots  were  taken  from  the  data 
points  of  Gunn  and  East  [1954]. 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


Ratio  of  reflectivities,  tjcl/’Ica 


^ 
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Figure  28.  Ratio  of  cloud  reflectivity  to  clear  air  irregularity 

reflectivity  as  a function  of  structure  constant,  gnd 
cloud  liquid  water  content  for  f = 15.7  GHz. 


Normalized  Attenuation  Coefficient,  k/M  (dB/km)/(gm/m3) 


Figure  29. 


t 

I 


Normalized  attenuation  coefficient  k /M^  for  water 
cloud  as  a function  of  frequency  and  te^erature. 
(Data  points  taken  from  Gunn  and  East  [1954].) 
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In  Fig.  30  are  plotted  further  theoretical  values  of  K for 
water  cloud  up  to  180  GHz*.  These  were  calculated  assuming  the  cloud 
drop  radii  are  very  much  smaller  than  the  wavelength  which,  in  general, 
is  the  case  over  the  frequency  range  considered. 

5.0  MULTIPLE  WAVELENGTH  RADAR  SYSTEMS  FOR  ESTABLISHING  PATH  ATTENUATION 
THROUGH  RAIN  AT  VARIABLE  FREQUENCIES 

A method  is  developed  here  for  using  a dual  wavelength  radar 
system  for  establishing  drop  size  distribution  aloft  from  which  necessary 
regression  relationships  at  arbitrary  wavelengths  above  10  GHz  may  be 
deduced.  These  regression  relations  may  subsequently  be  used  to  establish 
path  attenuations  using  an  S-band  radar.  This  multiple  wavelength  system 
represents  a modification  of  a previously  examined  dual  wavelength  method 
(Goldhirsh  and  Katz  [1974])  and  removes  some  of  the  ambiguities  arising 
in  this  former  system. 

5. 1 Rain  Drop  Size  Distribution  Using  a Dual  Wavelength  Method 

A technique  has  been  suggested  by  several  investigators  for 
using  a dual  wavelength  system  to  ascertain  the  drop  size  distribution 
aloft  using  multiple  wavelength  radars  (Goldhirsh  and  Katz  [1974],  Atlas 
and  Ulbrich  [1974]).  Using  the  measured  drop  size  distributions  associ- 
ated with  the  pulse  volumes  along  a slant  path,  the  individual  attenuation 
coefficients  at  arbitrary  frequencies  may  be  conq^uted  and  the  overall 
attenuation  established.  Goldhirsh  and  Katz  [1974]  and  Goldhirsh  [1975a] 
have  made  error  analyses  for  the  wavelength  combinations  (1  cm  - 10  cm), 

(3  cm  - 10  cm),  and  (1  cm  - 3 cm)  and  found  the  (X  = 1 cm,  A = 10  cm) 
wavelength  combination  least  sensitive  to  measurement  uncertainties 
altliough  range  limitation  problems  may  arise  due  to  severe  attenuation 
for  the  A = 1 cm  signals  at  the  heavier  rainfalls. 

A distinct  advantage  of  the  dual  wavelength  technique  is  that 
no  empirical  relation  of  the  form  aZ*’  needs  to  be  assumed  a priori.  A 

* Results  are  contained  in  CCIR  Study  Group  Doc  2/181E,  Doc.  5/161-E, 

9 May  1977,  Draft  Report  234-3  (Rev.  76). 


, (dB/km)/(gm/m'^) 
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Figure  30.  Theoretical  values  of  K=  /M<.  for  water  cloud  up  to 
180  GHz.  These  results  are  contained  in  CCIR  Study 
Group  Doc.  2/181-E,  Doc.  5/161-E,  9 May  1977,  Draft 
Report  234-3  (Rev.  76). 
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basic  assumption  in  this  method  is  that  the  drop  size  distribution 
possesses  the  general  exponential  form, 

N(D)dD  = N^exp(-AD)dD 


(5.1) 


where 

N(D)dD  = the  number  of  drops  of  diameters  between  D and  D+dD 
per  unit  volume 

N^,A  = variable  parameters  of  the  exponential  form  (5.1) 
determined  by  the  dual  wavelength  method 

The  validity  of  this  assumption  has  been  examined  by  Goldhirsh  [1976c]  and 
Ulbrich  and  Atlas  [1977]  and  found  using  their  samples  sizes  to  be  very 
good. 

By  way  of  review,  the  method  assumes  power  returns  from  two 

coincident  pulse  volumes;  one  from  a radar  operating  at  a non- attenuating 

wavelength,  P^(r),  and  the  other  from  a radar  operating  at  an  attenxiating 

wavelength,  P^(r),  where  r is  the  range  to  the  pulse  volume  center  (Fig. 

31).  Hence,  P (r)  and  P (r)  are  given  by, 

O 3i 


p 

o 


(r) 


Co  Z(r) 
2 


r 

C n(r)  -0.2  f k dr 

P^(r)  = 10  o*' 

r 


(5.2) 


(5.3) 


where  C^  and  are  the  appropriate  radar  constants,  Z is  the  reflectivity 
factor,  and  n is  the  Mie  reflectivity  at  the  attenuating  wavelength.  We 
express  the  power  return  from  the  range,  r ♦ A,  in  the  same  form  as  (5.3) 
•nii  jerive  by  straightforward  methods  the  relation, 


2 


k dr 


(5.4) 
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from  which. 


k = 


P^Cr)  n(r+A) 
(r+A)^  P^(r+A)  n (r) 


where  A [km]  is  the  range  averaging  interval  (Fig.  31)  and  k .1  s the 
average  attenuation  coefficient  over  the  interval,  A. 

An  alternate  form  of  (5.5)  may  be  derived  by  dividing  (5.3)  by 
(5.2)  at  the  different  range  intervals.  We  arrive  at. 


k = 


Pg^(r)  PQ(r+A)  n(r+A)  Z(r) 
P^(r+A)  P^(r)  n(r)  Z(r+A) 


which  is  range  normalized.  We  note  that  k is  a directly  measurable 
quantity  expressed  simply  in  terms  of  the  ratio  of  the  powers  only  when 
either  of  the  following  conditions  hold, 

n(r)  = n(r+A) 


or 


ri(r-t-A)  „ Z(r-t-A) 
n(r)  - Z(r) 

The  equality  (5.7)  implies  the  rain  is  uniform  at  the  range  bins  centered 
at  r and  r+A,  although  it  may  be  non-uniform  in  the  intervening  bins.  An 
examination  of  the  validity  of  this  equality  at  various  wavelengths  and 
rain  rate  combinations  will  be  examined  shortly. 

Considering  (5.2),  the  average  reflectivity  factor,  Z,  over  the 
range  interval.  A,  is  given  by, 

r+A 

2 

o J 

r 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


(5.9) 
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The  reflectivity  factor,  Z,  and  attenuation  coefficient,  k,  may 
be  expressed  theoretically  by. 


u 

-"o  / 


max 


D exp(-AD)dD 


min 


u 

/ 


min 


max 


ext 


(D)  exp(-AD)dD 


where  the  exponential  drop  size  distribution  form  given  by  (5.1)  is 

assumed  in  the  general  forms,  (2.28)  and  (2.29).  In  the  above  equations, 

D^.  , D are  minimum  and  maximum  drop  size  diameters  of  the  distribu- 
min  max  ‘ 

tion  (e.g.,  .03  and  0.5  cm,  respectively). 

If  the  conditions  (5.7)  or  (5.8)  are  valid,  Z and  k in  (5.10) 
and  (5.11)  are  radar  measurable  quantities  and  those  equations  have  the 
unknowns  and  A which  may  be  solved  numerically.  Hence  an  equivalent 
exponential  distribution  having  associated  with  it  a unique  pair  of 
values,  (N^,  A)  may  be  established  for  the  interval,  A.  The  attenuation 
coefficient,  k^,  at  any  other  wavelength  may  be  calculated  from  (5.11) 
for  each  interval,  A^,  and  the  integrated  one  way  attenuation  obtained 
using. 


i=l 


where  n represents  the  total  number  of  path  intervals,  A. 

5.1.1  Errors  Due  to  Inhomogeneous  Rain  Rates 

We  examine  here  the  error  in  IT  due  to  inhomogeneous  rain  rates 
along  the  path.  This  error  arises  when  (5.8)  is  assumed  in  (5.6)  and  k 


(5.10) 


(5.11) 


(5.12) 
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is  expressed  simply  in  terms  of  a ratio  of  received  powers.  In  Fig.  32, 
the  Mie  reflectivity  [m  is  plotted  as  a function  of  rain  rate  for 

X-band  (10  GHz),  K -band  (15.7  GHz),  and  K -band  (35  GHz).  These  curves 

u di 

were  deduced  from  Haddock  [1956]  who  assumed  a Laws  and  Parsons  distri- 
bution. Also  plotted  is  the  Rayleigh  reflectivity  at  3 GHz  as  a function 
of  rain  rate  arrived  at  using  (2.9)  and  Z = 200  R^'^. 

Using  the  reflectivity  results  in  Fig.  32,  the  attenuation 
coefficient  error  defined  as. 


n(r-t-A)  Z(r) 
n(r)  Z(r+A) 


is  tabulated  in  the  form  2A  6k  in  Table  5.1  for  various  rain  rate  combi- 
nations at  r and  r+&.  Also  tabulated  are  the  values  of  k averaged  over 
the  indicated  rain  rates  for  the  various  frequencies.  In  arriving  at 
these  averages,  the  attenuation  coefficient  versus  rain  rate  relations 
were  obtained  by  interpolating  the  tabulations  of  Medhurst  [1965].  These 
results  are  plotted  in  Fig.  33.  In  the  averaging  process,  an  exponential 
variation  of  rain  rate  with  range  over  the  interval.  A,  was  assumed.  In 
Table  5.2  are  tabulated  the  corresponding  percentage  errors,  6k/k  x 100, 
for  the  various  rain  rate  intervals  and  frequency  combinations.  In  these 
tabulations  it  is  assumed  2A  = 1 km.  For  other  intervals  A,  the  indicated 
error  values  should  be  divided  by  2A.  For  example,  for  A = 0.15  km,  the 
errors  in  estimating  k (from  Table  5.2)  in  the  rain  rate  interval  40  to 
SO  mm/hr  is  97,9"o,  56. 9»,  and  51.6%  for  the  X-S,  K -S,  and  K -S  band 
cases,  respectively  (i.e.,  3.33  times  the  corresponding  values  in 

Table  5.2).  We  note  that  although  the  attenuation  errors  6k  are  smallest 
at  X-band,  the  percentage  errors  are  largest. 

The  above  results  hence  give  sound  notions  as  to  the  magnitude 
of  errors  that  may  be  introduced  in  k for  the  various  wavelength  combi- 
nations due  to  the  assumption  (5.8).  These  may  translate  into  sizeable 
errors  in  the  estimation  of  drop  size  distribution  and  ultimately  in 
attenuation  coefficient  at  other  wavelengths  (see,  for  example,  Goldhirsh 
and  Katz  [1974]  and  Goldhirsh  [1975]). 


Reflectivity,  tj  (S-band),  m 


f=  15.7  GHz  f = 3GHz 
(Ku  band)  (S-band) 


f = 10GHz 
(X-band) ^ 


f = 35GHz 


Reflectivity,  t\  (X,  K^,  Ka  bands),  nr 


Figure  32.  Mie  reflectivity  as  a function  of  rain  rate  for  X-band 
(10  aiz),  Ku-band  (15.7  Qlz),  and  Kg-band  (35  GHz) 
assuming  a Laws  and  Parsons  distribution.  Also  plotted 
is  t)ie  S-band  (3  GHz)  Rayleigh  reflectivity. 
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Table  5.1 

TABLE  OF  VALUES  OF  k AVERAGED  OVER  THE  INDICATED  RAIN  RATE 
INTERVALS  AND  ASSOCIATED  REFLECTIVITY  RATIO  ERRORS,  2A  6k. 
(X  = 10  GHz,  u = 15.7  GHz,  a = 35  GHz) 


1 R(r) 

R(r+A) 

k 

X 

k 

u 

k 

a 

2A  6k  (dB) 

iran/hr 

dB/km 

X-S 

K -S 
u 

1 K -S 

1 a 

2.5 

5.0 

0.0 

0.9 

0.0 

0.0 

B 

5.0 

10.0 

0.1 

0.0 

0.1 

mm 

10.0 

20.0 

0.3 

0.0 

0.5 

BB 

20.0 

40.0 

0.7 

0.3 

0.6 

40.0 

80.0 

mgm 

12.9 

0.5 

0.7 

80.0 

160.0 

B 

9.0 

25.4 

0.8 

1.25 

B 

Table  5.2 

TABLE  OF  PERCENTAGE  ERRORS  DUE  TO  REFLECTIVITY 
RATIOS  OVER  INDICATED  ^ 

RAIN  RATE  INTERVAL  ON  A PER  km  BASIS  (i.e.,  ~ » 100). 


R(r) 

R(r+A) 

(6k/k)  X 100  (%) 

mm/hr 

X-S 

(10-3  GHz) 

K -S 
u 

(15.7-3  GHz) 

V" 

(35-3  GHz) 

2.5 

5 

-- 

111.1 

5 

10 

-- 

25.0 

94.1 

10 

20 

-- 

55.5 

51.5 

20 

40 

42.8 

31.6 

27.3 

40 

80 

29.4 

17.1 

15.5 

80 

160 

19.0 

13.9 

9.5 

R (mm/hr) 
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It  may  be  noted  that  the  errors  due  to  the  assumption  (5.7)  in 
presence  of  inhomogeneous  rain  rates  have  not  been  considered  here  as 
they  are  more  severe  than  those  arising  from  the  assumption  (5.8). 

In  the  next  paragraphs  are  described  a modified  dual  wavelength 
method  which  may  ultimately  reduce  the  above  type  errors  and  be  more 
practical  for  implementation. 

5.2  Modified  Dual  Wavelength  Method 

We  proceed  here  on  the  basis  that  a dominant  drop  size  distri- 
bution may  be  assumed  over  at  least  each  rain  period  and  that  the  corre- 
sponding statistics  can  be  deduced  using  (a)  the  dual  wavelength  method 
applied  only  over  those  intervals.  A,  having  uniform  rain  rates;  that  is, 
where  the  equality  (5.7)  is  valid  and/or  (b)  a ground  disdrometer  as,  for 
example,  the  Rowland  [1975,  1976]  or  Joss  [1976]  types. 

For  each  sampling  of  the  drop  size  distribution,  N(D),  a value 
of  the  Rayleigh  reflectivity  factor,  Z,  attenuation  coefficient,  k,  at 
the  attenuating  wavelength  of  the  radar,  and  the  corresponding  Mie 
reflectivity,  n,  may  be  determined  from  the  equations. 


•max 

Z = # N(D)dD  [mmVm'^] 


/■ 

min 


(5.14) 


/max 

^ext^'^^  N(D)dD  [dB/km] 

u . 
min 


(5.15) 


u 

/■ 

D 

min 


a(D)  N(D)dD 


[m‘^] 


where  all  of  the  above  parameters  have  been  previously  defined  except 

2 

o(D)  which  represents  the  Mie  scattering  cross  section  [m  ] and  is 


(5.16) 
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dependent  upon  drop  diameter,  wavelength,  and  refractive  index.  This 
quantity  has,  for  example,  been  tabulated  by  Stephens  [1961]  as  a 
function  of  drop  diameter  for  a series  of  wavelengths. 

With  k,  Z,  and  n calculated  for  each  drop  size  distribution, 
k-Z  and  n-Z  scatter  diagrams  may  be  arrived  at  from  which  best  fit 
regression  relations  of  the  form, 

k = a 2 

^ -7B 

n = A z 

The  next  step  in  the  modified  dual  wavelength  method  is  to 
check  (5.17)  and  (5.18)  against  the  attenuating  wavelength  measurement 
(5.4)  in  all  regions  of  rain  including  those  intervals.  A,  where  the  rain 
is  non-uniform.  This  may  be  referred  to  as  the  drop  size  distribution 
"truth"  measurement  phase  of  the  method.  Hence  substituting  (5.17)  and 

(5.18)  into  (5.4)  we  obtain  for  the  ratio  of  powers  at  the  attenuating 
wavelength. 


10  Log^Q 


r P (r) 
a 


(r+A)  P (r+A) 
& 


= 10  B Log 


10 


lilL 


Z(r+A) 


+ 2a 


r+A 

/ 


Z dr 


The  right  hand  side  (RHS)  of  (5.19)  may  be  obtained  from  the  S-band 
measurement  of  Z(r)  and  a straightforward  calculation.  The  left  hand 
side  (LHS)  is  arrived  at  through  a relative  power  measurement  at  the 
attenuating  wavelength.  Comparison  of  the  LHS  and  RHS  of  (5.19)  in  both 
range  and  time  enables  one  to  check  the  validity  of  using  (5.17)  and 

(5.18)  obtained  from  the  original  drop  size  distribution.  It  is,  of 
course,  not  reasonable  to  expect  a one-to-one  check  for  all  ranges  and 
times.  However,  if  a dominant  drop  size  distribution  does  exist,  then 
the  check  should  at  least  show  statistically  good  results.  The  LHS  of 

(5.19)  represents  a relative  measurement  and  therefore  could  be  acquired 
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with  good  accuracy  (e.g.,  0.5  dB  uncertainty).  The  RHS  could  be  in  error 
due  to  the  variances  associated  with  a,  b.  A,  and  B or  due  to  the  uncer- 
tainty in  the  measurement,  Z,  which  requires  an  absolute  calibration.  In 
the  event  the  RllS  systematically  underestimates  or  overestimates  the  UlS, 
the  DSD  may  be  adjusted  such  that  the  RllS  produces  a better  check  for  the 
revised  enq^irical  parameters.  A technique  for  taking  out  the  uncertainty 
in  k is  described  in  Section  5.3.3. 

Once  confidence  in  the  proper  distribution  is  established 
through  repeated  checks  using  (5.19),  only  the  S-band  radar  need  be  used 
and  the  reflectivity  profile,  Z(r),  established.  The  attenuation  at  any 
other  wavelength  is  merely. 


I 

A = ^ a Z*^(r)  dr  [dB] 

o*^ 

where  a and  b at  other  wavelengths  are  derived  through  scattergrams  of 
k versus  Z computed  using  the  original  (or  corrected)  drop  size  distribu- 
tion data  base. 

In  summary,  the  steps  involved  in  using  the  modified  diial  wave- 
length method  are : 

(1)  Measure  the  power  scattered  from  rain  at  both  S-band 

(3  GHz)  and  either  K^-band  (15.7  GHz)  or  K^-band  (35  GHz) 
(the  latter  representing  a preferable  wavelength)  using 
coincident  beams. 

(2)  Over  range  intervals.  A,  only  where  the  reflectivities  are 
uniform,  the  dual  wavelength  method  should  be  implemented 
and  samples  of  N(D)  acquired.  If  possible,  ground  measure- 
ments of  DSD  should  also  be  obtained  with  disdrometers. 

(3)  Using  the  measured  S2unples  of  DSD,  Z,  k,  and  n as  given  by 
(5.14),  (5.15),  and  (5.16),  respectively,  should  be  cal- 
culated; k and  n correspond  to  the  attenuating  wavelength 
of  the  radar. 


(5.20) 
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(4)  From  the  results  of  (3),  scatter  diagrams  of  log  k and 
log  n versus  log  Z may  be  plotted  and  best  fit  empirical 
relations  (5.17)  and  (5.18)  at  the  attenuating  wavelength 
of  the  radar  may  be  determined. 

(5)  The  empirical  parameters,  a,  b,  and  B and  the  S-band 
reflectivity  levels  Z(r)  may  now  be  injected  in  the  RHS 
of  (5.19)  and  checked  against  the  relative  power  measure- 
ments of  the  attenuating  wavelength  radar  [LHS  of  (5.19)] 
for  all  rain  rate  ranges. 

(6)  If  the  LHS  under-  or  overestimates  the  RHS,  equivalence 
of  the  two  sides  may  be  achieved  by  methods  described  in 
Section  5.3.3. 

(7)  With  the  RHS  and  LHS  showing  good  statistical  checks  for  a 
fixed  calibration  constant  of  the  S-band  radar,  a dominant 
DSD  is  therefore  established  from  which  scatter  diagrams 
of  Log  k versus  Log  Z may  be  obtained  at  other  attenuating 
wavelengths  where  knowledge  of  the  path  attenuation  is 
desired. 

(8)  Forming  empirical  relations  of  the  form  (5,17)  from  (7), 
appropriate  parameters  a and  b may  be  determined  at  other 
wavelengths. 

(9)  The  empirical  parameters  a and  b may  then  be  used  in  con- 
junction with  the  S-band  reflectivity  factor  measurements, 
Z,  to  obtain  the  integrated  path  attenuation  as  given  by 
(5.20)  for  other  wavelengths. 

5.2.1  Pertinent  Parameters  at  15.7  and  35  GHz. 

In  this  section  are  provided  the  necessary  data  to  implement 
the  modified  dual  wavelength  method  at  either  the  frequencies  of  15.7  or 
35  GHz.  In  Fig.  34  are  plotted  k/N^  versus  A at  these  frequencies 
obtained  through  the  numerical  integration  of  the  integral. 
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U 

f 

<<0  J 


max 


min 


where 


Cext(D)  exp(-AD)dD 


[dB/km] cm 


D = 0 . 5 cm 
max 


D . = 0.03  cm 

min 


and  the  values  of  are  obtained  by  interpolating  the  tabulations  in 

Medhurst  [1965J.  The  resulting  tabulations  of  vs  D are  given  in 

Table  5.3.  We  note  that  the  units  of  N are  cm”'*. 

o 

The  indicated  minimum  and  maximum  diameters  in  (5.22)  were 
used  as  these  were  found  to  bound  the  drop  diameters  of  most  rain  rates 
(Goldhirsh  [1976c]).  Also  plotted  in  Fig.  34  is  the  normalized  rain  rate 

4 

R/N^  [(mm/hr) cm  ] as  a function  of  A.  This  curve  was  obtained  through 
the  numerical  integration  of, 


^ = 1.885  X 10^ 


u 

J 


min 


max 


v(D)  D exp(-AD)dD 


[mm/hr] (cm  ) 


where  v(D)  is  the  drop  terminal  velocity  (m/sec)  and  D is  expressed  in 
cm.  The  tabulated  data  given  in  Table  5.4  was  used  in  the  numerical 
integration  of  (5.23). 

In  Fig.  35  are  plotted  the  normalized  Mie  reflectivity  n/N 

-14  ° 

(m  ) (cm  ) at  15.7  and  35  GHz.  These  curves  were  obtained  through  the 

numerical  integration  of. 


^ = 10‘ 


L 

I 


min 


max 


a(D)  exp(-AD)dD 


(m‘^(cm^) 


(5.21) 


(5.22) 


(5.23) 


(5.24) 
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Table  5.3 

TABUUTIONS  OF  Cg^t  [(dB/km)cm3]  vs  D [cm]  OBTAINED 

BY  INTERPOLATING  THE  RESULTS  OF  MEDHURST  [1965]. 

■ 

D 

(dB/km)cm^ 

cm 

f ■ 35  GHz 

0.03 

8.90  + 0 

1.89  + 1 

0.05 

1.05  + 1 

6.65  ♦ 1 

0.10 

1.60  + 2 

1.54  + 3 

0.15 

1.33  + 3 

9.03  + 3 

0.20 

6.83  + 3 

2.49  + 4 

0.25 

1.99  + 4 

5.06  + 4 

0.30 

3.23  + 4 

8.61  + 4 

0.35 

5.25  + 4 

1.23  + 5 

0.40 

8.11  + 4 

1.58  + 5 

0.45 

1.21  + 5 

1.86  + 5 

0.50 

1.71  + 5 

2.36  5 

Table  5.4 

TABULATIONS  OF  THE  INTEGRAND  FACTOR 
OF  RAIN  RATE  TAKEN  FROM  MEDHURST  [1965], 


D 

cm 

V 

ra/sec 

< •«  ' 

1.885  X 10^  V 

0.03 

... 

1.60  + 2 

0.05 

2.06 

4.85  + 2 

0.10 

4.03 

7.60  + 3 

0.15 

5.40 

3.44+4 

0.20 

6.49 

9.79  + 4 

0.25 

7.41 

2.18  + 5 

0.30 

8.06 

4.10  + 5 

0.35 

8.53 

6.89+5 

0.40 

8.83 

1.07  + 6 

0.45 

9.00 

1.55  + 6 

0.50 

i 

9.09 

2.14  + 6 

106  _ 
E 

u 


{f  = 35GH2) 
No 


^ If  = 15.7  GHz) 
Nq 


o 
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Figure  35.  Normalized  Mie  reflectivities  at  15.7  and  35  GHz  plotted 
as  a function  of  A.  Normalized  reflectivity  factor,  Z, 
as  a fmction  of  A also  shown. 
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2 

where  o(D)  [cm  ] is  the  Mie  scattering  cross  section  of  a drop  of 
diameter,  D [cnj]  for  the  given  wavelength.  The  scattering  cross  sections 
were  obtained  from  the  tabulation  of  Stephens  [1961]  with  the  values  at 
15.7  interpolated.  These  are  listed  in  Table  5.5. 

Also  plotted  in  Fig.  35  is  the  normalized  reflectivity  factor 
Z 6 3 4 

jj—  (mm  /m  )cm  versus  A.  This  was  obtained  through  the  integration  of, 

‘^o 

*^max 

|-=(lxlO^^)  exp(-AD)dD  [ (mmW)  (cm"^)  ] 

min 


It  is  interesting  to  note  that  the  righthand  side  of  (5.25)  may  be 
analytically  integrated  and  is  given  by. 


Z 

N 


f(A, 


D ) 

max"^ 


f(A,  D . ) 
^ ’ min-’ 


where , 


6 

f(A,  D)  s 10^^  exp(-AD) 

p=0 


(-l)P 


61  D^-P 
(6-p)!  (.A)P*^ 


In  Table  5.6  are  tabulated  the  numerically  calculated  values  of 

R/N^,  k/N^,  and  n/N^  as  a function  of  A plotted  in  Figs.  34  and  35.  In 

Table  5.7  are  tabulated  Z/N  versus  A for  the  cases  in  which  D = 0.5 

0 max 

cm  and  ■».  We  note  that  this  latter  quantity  (i.e..  Z/N  for  D = oo^ 

^ ' o max  ■' 

which  is  often  used  for  sinplicity  overestimates  the  actual  reflectivity 
factor.  Also  tabulated  for  reference  is  the  Marshall-Palmer  [1948]  rain 
rate  as  a function  of  A given  by, 

/4n4.76 

R = [mm/hr] 

where  A is  in  cm"^. 


L 


(5.25) 


(5.26) 


(5.27) 


(5.28) 
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Table  5.5 

TABULATIONS  OF  MIE  SCATTERING  CROSS  SECTION,  a [cm^]  VERSUS 
RAINDROP  DIAMETER,  D,  OBTAINED  FROM  STEPHENS  [1961]. 


Table  5.6 


TABULATIONS  OF  R/Nq,  K/Nq,  AND  n/N©  VERSUS  A OBTAINED 
FROM  THE  NUMERICAL  INTEGRATIONS  OF  (5.23),  (5.21),  AND  (S,24),  RESPECTIVELY. 


A 

m)‘^ 

R/N 

o 

z 

, (dB/km) (cm) 

4 

K 

1 

■ 

4 

f = 15. 

7 

GHz 

f = 35  GHz 

BBj 

mm 

f = 35 

GHz 

1 

1.74 

+ 

5 

1.34 

+ 

4 

2.61 

4 

1.68 

2.22 

2 

1.16 

4- 

5 

9.01 

+ 

5 

1.79 

+ 

4 

1.11 

1.64 

3 

7.89 

4 

6.09 

+ 

3 

1.24 

+ 

4 

7.31 

1 

1.22 

4 

5.38 

+ 

4 

4.14 

+ 

3 

8.66 

+ 

4 

4.85 

1 

9.06 

- 

1 

5 

3.70 

•f 

4 

2.84 

+ 

3 

6.10 

+ 

3 

3.23 

1 

6.78 

- 

1 

6 

2.56 

4- 

4 

1.95 

+ 

3 

4.34 

3 

2.16 

1 

5.10 

- 

1 

7 

1.79 

+ 

4 

1.36 

+ 

3 

3.11 

+ 

3 

1.45 

1 

3.85 

- 

1 

8 

1.26 

4* 

4 

9.53 

+ 

2 

2.25 

+ 

3 

9.85 

2 

2.93 

- 

1 

9 

9.00 

4- 

3 

6.73 

+ 

2 

1.65 

+ 

3 

6.71 

2 

2.24 

- 

1 

10 

6.49 

4* 

3 

4.79 

+ 

2 

1.21 

+ 

3 

4.60 

2 

1.71 

- 

1 

11 

4.71 

+ 

3 

3.44 

+ 

2 

9.04 

+ 

2 

3. 16 

2 

1.31 

- 

1 

12 

5.46 

4* 

3 

2.49 

+ 

2 

6.79 

+ 

2 

2.20 

2 

1.02 

- 

1 

13 

2.58 

4- 

3 

1.83 

+ 

2 

5.16 

+ 

2 

1.54 

2 

7.89 

- 

2 

14 

1.94 

4- 

5 

1.34 

+ 

2 

3.95 

+ 

2 

1.08 

2 

6.15 

- 

2 

15 

1.46 

4- 

3 

9.94 

+ 

1 

3.05 

2 

7.63 

• 

3 

4.81 

- 

2 

16 

1.13 

4* 

3 

7.44 

+ 

1 

2.39 

+ 

2 

5.44 

w 

3 

3.78 

- 

2 

17 

8.71 

4- 

2 

5.60 

+ 

1 

1.88 

+ 

2 

3.90 

m 

3 

2.98 

- 

2 

18 

6.81 

4- 

2 

4.26 

+ 

1 

1.49 

+ 

2 

2.81 

3 

2.36 

- 

2 

19 

5.39 

4- 

2 

3.26 

+ 

1 

1.19 

+ 

2 

2.05 

3 

1.88 

- 

2 

20 

4.29 

4- 

2 

2.51 

+ 

1 

9.54 

+ 

1 

1.50 

3 

1.49 

- 

2 

21 

3.45 

4- 

2 

1.95 

+ 

1 

7.73 

+ 

1 

1.11 

3 

1.20 

- 

2 

22 

2.80 

4* 

2 

1.51 

+ 

1 

6.30 

+ 

1 

8.21 

■ 

9.60 

- 

3 

23 

2.29 

4- 

2 

1.19 

+ 

1 

5.16 

+ 

1 

6. 14 

■ 

7.73 

- 

3 

24 

1.89 

4- 

2 

9.39 

4.26 

1 

4.61 

■ 

6.25 

- 

3 

25 

1.56 

4- 

2 

7.45 

3.54 

♦ 

1 

3.49 

■ 

5.08 

- 

3 

26 

1.30 

4“ 

2 

5.94 

2.95 

+ 

1 

2.66 

■ 

4.13 

- 

3 

27 

1.09 

♦ 

2 

4.76 

2.48 

1 

2.04 

■ 

3.36 

- 

3 

28 

9.23 

4* 

1 

3.84 

2.09 

+ 

1 

1.58 

■ 

2.76 

- 

3 

29 

7.83 

4- 

1 

3.10 

1.76 

1 

1.22 

■ 

2.26 

- 

3 

30 

6.69 

4- 

1 

2.53 

i 

1.50 

1 

9.54 

m 

5 

1.86 

- 

3 

31 

5.73 

4- 

1 

2.06 

1.29 

♦ 

1 

7.49 

m 

5 

1.54 

3 

32 

4.94 

+ 

1 

1.69 

1.10 

1 

5.93 

m 

5 

1.28 

- 

3 

33 

4.28 

4- 

1 

1.40 

9.45 

4.70 

m 

5 

1.06 

- 

3 

34 

3.71 

4- 

1 

1.16 

8.16 

3.76 

m 

5 

8.86 

- 

4 

35 

3.24 

♦ 

1 

9.60 

- 

1 

7.06 

3.03 

m 

5 

7.41 

- 

4 

36 

2.84 

4- 

1 

8.05 

- 

1 

6.14 

2.44 

m 

5 

6.21 

- 

4 

37 

2.49 

4- 

1 

6.73 

- 

1 

5.35 

1.99 

m 

5 

5.23 

- 

4 

38 

2.19 

4 

1 

5.66 

- 

1 

4.69 

1.63 

m 

5 

4.40 

- 

4 

39 

1.94 

4 

1 

4.80 

- 

1 

4.10 

1.34 

m 

5 

3.71 

- 

4 

40 

1.73 

4 

1 

4.08 

- 

1 

3.61 

1.10 

m 

s 

3.15 

- 

4 
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Table  5.7 

TABULATIONS  OF  Z/N^  FOR  = 0.5  an  AND  = «. 
ALSO  PLOTTED  ARE  THE  MARSHALL- PAIAIER 
RAIN  RATES  [1948]  AS  A FUNCTION  OF  A. 


B 

Z/N^  - (nun^/m^)  cm^ 

R (M-P) 
(inm)/hr 

D = 0.5  on 

max 

D S 00 

max 

7.24  + 8 

7.20  + 14 

m m 

2 

4.68  + 8 

5.62  + 12 

3 

3.05  + 8 

3.29  + 11 

4 

1.99  + 8 

4.39  + 10 

5 

1.31  + 8 

9.22  + 9 

6 

8.62+7 

2.57  + 9 

« ^ 

7 

5.71  + 7 

8.74  + 8 

8 

3.79  + 7 

3.43  + 8 

9 

2.54  + 7 

1.50  + 8 

10 

1.71  + 7 

7.20+7 

8.26  + 2 

11 

1.16  + 7 

3.69  + 7 

5.25  + 2 

12 

7.83  + 6 

2.01  + 7 

3.47  + 2 

13 

5.43  + 6 

1.15  + 7 

2.37  + 2 

14 

3.76  + 6 

6.83  + 6 

1.66  + 2 

15 

2.62  + 6 

4.21  + 6 

1.20  + 2 

16 

1.84  + 6 

2.68  + 6 

8.82  + 1 

17 

1.31  + 6 

1.75  + 6 

6.61  + 1 

18 

9.33  + 5 

1.18  + 6 

5.03  + 1 

19 

6.75  + 5 

8.05+5 

3.89  + 1 

20 

4.89+5 

5.63  + 5 

3.05  + 1 

21 

3.59  + 5 

4.00+5 

2.42  + 1 

22 

2.66  + 5 

2.89+5 

1.94  + 1 

23 

1.99  + 5 

2.12  + 5 

1.57  + 1 

24 

1.50  + 5 

1.57  + 5 

1.28  + 1 

25 

1.14  + 5 

1.18  + 5 

1.05  + 1 

26 

8.73+4 

8.96  + 4 

8.7 

27 

6.75  + 4 

6.88  + 4 

7.3 

28 

5.26  + 4 

5.34  + 4 

6.1 

29 

4.13  + 4 

4.17  + 4 

5.2 

30 

3.27  + 4 

3.29+4 

4.4 

31 

2.60  + 4 

2.62  + 4 

3.8 

32 

2.09+4 

2.10  + 4 

3.3 

33 

1.68  + 4 

1.69  + 4 

2.8 

34 

1.37  + 4 

1.37  + 4 

2.4 

35 

1.12  + 4 

1.12  + 4 

2.1 

36 

9.18  + 3 

9.18  + 3 

1.9 

37 

7.58  + 3 

7.58  + 3 

1.6 

38 

6.29  + 3 

6.29  + 3 

1.4 

39 

5.24  + 3 

5.24  + 3 

1.3 

40 

4.39+3 

4.39  + 3 

1.1 
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5.2.2  Solving  for  the  Drop  Size  Distribution 

In  solving  for  and  A using  the  dual  wavelength  method,  it  is 
convenient  to  use  the  ratio. 


k 

Z 


u 

/ 


min 


max 


Cg^^(D)  exp(-AD)dD 


f(A,D  ) - f(A,D  . ) 
'■  ’ max-^  '■  ’ min-^ 


(5.29) 


since  it  is  independent  of  and  is  only  a function  of  A-  Hence  a 

measurement  of  k and  Z uniquely  defines  a value  of  A,  using  (5.29).  The 

above  ratios  may  be  obtained  from  the  tabulations  of  k/N^  and  Z/N^  in 

Tables  5.6  and  5.7  and  are  plotted  in  Figure  36  for  f = 15.7  and  35  GHz 

as  a function  of  A.  With  A known,  the  appropriate  value  of  k/N  or  Z/N^ 

may  be  obtained  (by  exponential  interpolation)  from  Tables  5,6  or  5.7  or 

from  Figures  34  and  35.  The  corresponding  value  of  is  calculated  by 

dividing  the  normalized  values  into  the  respective  measured  values  (k  or 

Z).  Hence,  the  pair  k and  Z enable  the  determination  of  N and 

o 

5.3  Error  Analysis  at  and  Bands 

We  examine  here  uncertainties  that  may  arise  in  the  modified 
dual  wavelength  method  due  to  realistic  uncertainties  in  the  measurement 
of  k and  Z from  which  empirical  parameters  are  obtained  using  the  dual 
wavelength  method.  Specifically,  we  assume  a ^ 2 dB  uncertainty  in  the 
measurement  of  Z and  a ^ 25%  error  in  k for  both  the  15.7  and  35  GHz 
cases.  With  these  uncertainties,  we  examine  the  degree  of  sensitivity  of 
the  RHS  of  (5.19)  to  these  errors  at  the  corresponding  frequencies.  It 
will  in  fact  be  demonstrated  that  the  uncertainties  at  35  GHz  is  far  more 
pronounced  that  at  15.7  GHz  and  that  this  enables  a correction  to  be  made 
in  both  the  calibration  of  the  attenuating  wavelength  measurements  and 
the  resulting  drop  size  distribution.  The  previously  used  measured  Z 
profile  from  which  empirical  rain  rates  were  deduced  in  Pig.  11  (solid 
curve)  is  used  here  as  a test  case.  The  partial  Z profile  measured  is 
given  in  Fig.  37. 
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Figure  36.  Ratio  of  attenuation  coefficient  to  reflectivity  factor 
as  a function  of  A for  f = 15.7  and  35  GHj. 


Measured  reflectivity  factor  Z (mm®/m3) 


Radar  range  (rm) 


Figure  37. 


Segment  of  reflectivity  profile  as  measured  at  Wallops 
Island,  Virginia,  on  July  10,  1976  at  19:52:20  GMT  at 
60*  azimuth  and  20°  elevation. 
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5.3.1  Comparison  of  the  Methods  at  15.7  and  35  GHz  Assuming 
Zero  Uncertainties 

To  arrive  at  a basis  for  comparison,  we  initially  examine  the 
RllS  of  (5.19)  assuming  a zero  uncertainty  using  the  Marshall-Palmer  value 
of  = 0.08  in  Tables  5.6  and  5.7.  The  resulting  values  of  k,  Z,  and  n 
enable  the  determination  of  best  fit  regression  relations, 

k = a (5.30) 

n = A Z®  (5.31) 


at  both  15.7  and  35  GHz.  The  corresponding  empirical  parameters  are 
given  in  Table  5.8. 

Substituting  the  appropriate  parameters  of  Table  5.8  into  the 
right  hand  side  of  (5.19)  and  using  the  example  Z profile  in  Fig.  37 
results  in  RHS  profiles  given  in  Figs.  38,  39,  and  40  for  averaging 
intervals  A = 0.15,  0.30,  and  0.60  km,  respectively.  We  note  that  the 
larger  interval  A = 0.6  km  improves  the  ability  to  distinguish  between 
the  two  frequencies  and  will  therefore  be  used  in  subsequent  error  com- 
parisons. We  also  note  in  Fig.  40  that  the  RHS  ranges  from  -9  to  +5  dB 
for  15.7  GHz  and  from  -4  to  +22  dB  for  35  GHz.  The  LHS  of  (5.19)  should 
be  measurable  to  within  ^0.5  dB  at  both  frequencies  since  it  represents 
a relative  measurement. 

5.3.2  Comparisons  of  the  Methods  at  15.7  and  35  GHz  Assuming 
Realistic  Measurement  Uncertainties 

We  here  assume  a ^ 2 dB  error  in  Z and  a ^ 25%  uncertainty  in 
the  measurement  of  k and  with  these  uncertainties  re-examine  the  RHS  of 
(5.19)  using  the  Z profile  of  Fig.  37  and  A = 0.6  km.  Hie  method  in 
which  this  error  analysis  is  performed  is  as  follows: 

(1)  k and  Z is  calculated  as  a function  of  A for  N =0.08 
o O 0 0 

using  Tables  5.6  and  5.7;  where  the  subscript  o in^lies 
the  "true  values". 


Table  5.8 

REGRESSION  PARAMETERS  a.  A,  b,  B OVER  INDICATED 
RAIN  RATE  INTERVALS,  R.  COEFFICIENT  OF  DETERMINATION, 
r2,  GIVES  MEASURE  OF  GOODNESS  OF  FIT. 


k(15.7) 

(dB/km) 

n(15.7) 

(m-1) 

k(35) 

(dB/km) 

n(35) 

(ro-1) 

a,  A 

3.25-4 

3.47  - 9 

5.48  - 3 

4.43  - 7 

b,  B 

0.835 

0.986 

0.685 

0.745 

2 

r 

0.99998 

0.99980 

0.99999 

0.99903 

R(inm/hr) 

10  - 206 

10  - 206 

1 

5 - 120 

5 - 120 

Range  normalized  power  ratio  IdB) 
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Figure  39.  Calculated  range  normalized  power  ratio  [RllS  of  (5.19)] 

for  f ■ 15.7  and  35  GHz  and  averaging  interval*  ^ = 0.3  km. 


Range  normalized  power  ratio  (dB) 


O f-  15.7  GHz 
□ f-35GHz 


8.8’ 9.0  9.2  9.4  9.6  aS  10.0  i6.2il6.410.6fl0'.8 11.011.2  1 1:4  11.6  1^ 

Radar  range  (km) 

Figure  40.  Calculated  range  normalized  power  ratio  [RUS  of  (5.19)]  for  f * 15.7 
and  35  GHz  and  averaging  interval,  A = 0.6  km.  Vertical  lines  depict 
extreme  uncertainties  assuming  a '•'25%  error  in  k and  ^2  dB  error  in  Z. 
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(2)  For  each  value  of  error  values  of  k and  Z are  calcu- 
lated. These  are  defined  as 

k+  = 1.25  k 

0 

k-  = 0.75  k 

0 

2+  = 1.58  Z 

0 

Z-  = 0.63  Z 

o 

(3)  Four  combinations  of  ratios,  k/Z  are  computed  for  each  A ; 
namely,  k+/Z+,  k-/Z+,  k+/Z-,  and  k-/Z-.  Using  the  curves 
of  Fig.  36,  error  values  of  A may  be  computed  for  each 
error  ratio. 

(4)  Using  the  four  error  values  for  A,  normalized  values  of 
k/N^  and  Z/N^  may  be  obtained  using  Tables  5.6  and  5.7  or 
Figs.  34  and  35. 

(5)  Four  error  values  of  may  be  obtained  by  dividing  the 

normalized  values,  and  Z/N^  into  the  respective  error 

values  of  k and  Z. 

(6)  Using  the  four  pairs  of  N^,  A calculated  corresponding  to 
each  error  combination,  we  may  calculate  four  error  values 
of  the  reflectivity,  n. 

D 

(7)  Best  fit  error  regression  relationships  n ^ A Z may  thus 
be  calculated  for  each  error  combination. 

(8)  Modified  best  fit  error  regression  relationships  of  k = 

b ^ 

Z^  ° may  also  be  calculated  for  each  error  combination. 

This  is  generally  given  by. 


k = 


b 

(1  - 6^)  ° 


a 


o 


(5.32) 


where  (uncertainties  in  k)  for  the  present  case  are 
^ 0.25  and  6^  (uncertainties  in  Z)  are  +0,68  and  -0.37. 


We  note  from  (5.32)  that  uncertainties  in  measured  values  of  Z 
and  k do  not  influence  b but  do  alter  a (to  be  demonstrated  shortly). 
Substituting  the  revised  values  of  a and  B into  the  RHS  of  (5.19),  new 
values  may  be  computed  for  each  error  combination. 

The  corresponding  error  values  of  the  RHS  as  a function  of 
range  for  each  error  combination  is  given  in  Table  5.9.  Also  given  are 
the  zero  error  cases  (subscript  o)  as  well  as  the  calculated  error  values 
of  the  empirical  parameters.  The  extreme  error  values  are  shown  plotted 
in  Fig.  40. 

It  is  interesting  to  note  from  Table  5.9  or  Fig.  40  that  the 
RHS  of  (5.19)  is  less  sensitive  to  the  given  errors  at  f »=  15.7  whereas 
at  f = 35,  significantly  greater  sensitivity  exists.  For  example,  we 
note  that  at  11.55  km  the  uncertainties  about  the  nominal  values  are  ^ 2 
dB  at  f = 15.7  Qlz  and  -6  to  +7  dB  at  35  GHz. 

The  larger  errors  in  the  K -S  band  case  in  Fig,  40  represent 

Si 

an  advantage  over  the  K^-S  band  case  in  that  these  can  be  used  to  reduce 
the  errors  in  the  original  measurement  in  k.  In  the  next  section  we  demon- 
strate how  an  analysis  of  the  differences  in  the  measured  left  hand  side 
and  the  calculated  right  hand  side  may  be  used  to  improve  the  measurement 
accuracy. 

In  utilizing  the  dual  wavelength  procedures  to  arrive  at  drop 
size  distributions  from  which  the  empirical  relation  k * a may  be 
deduced  at  other  wavelengths,  significantly  smaller  errors  may  be  shown 
to  arise  using  the  K^-S  band  case  as  compared  to  the  K^-S  band  configur- 
ation. To  demonstrate  this,  we  show  plotted  in  Fig.  41  worst  case  pre- 
dicted k-Z  curves  at  f » 15.7  GHz  predicted  using  the  dual  wavelength 
method  for  the  K -S  band  case  (35-3  GHz)  assuming  a +25%  error  in  the 
K^-band  measurement  and  a ^2  dB  error  in  the  Z measurement.  TTie  per- 
centage error  in  k for  Z values  below  Z = 10^  are  in  general  well  below 
^ 40%.  On  the  other  hand,  in  Fig.  42  are  shown  a similar  set  of  curves 
at  f = 35  GHz  predicted  from  the  K^-S  band  case  (15.7-3  GJlz).  The  upper 
error  case  results  in  errors  in  k of  more  than  a few  hundred  percent  and 


Attenuation  coefficient,  k (dB/km) 


Figure  41. 


Note:  k = a 2“ 

f = 15.7  GHz  predicted 
from  ka-S  band  system 


4 6 8 10^  2 4 6 8 10® 

Reflectivity  factor,  Z (mm®/m3) 


Best  fit  predicted  attenuation  coefficient  (k)  - reflectivity 
factor  (Z)  curves  for  f = 15.7  GHz  using  the  dual  wavelength 
method  for  the  band  case  (35-3  GHz).  Upper  and  lower 
error  curves  correspond  to  extreme  errors  assuming  *2  dB 
uncertainty  in  Z and  ^25%  uncertainty  in  k.  The  zero  error 
case  corresponds  to  Nn  ■ .08  cm*^. 


Atteunation  coefficient,  k (dB/km) 
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Reflectivity  factor,  Z (mm®/m3) 


Best  fit  predicted  attenuation  coefficient  (k)  - reflectivity 
factor  (Z)  curves  for  f * 35  GHz  using  the  dual  wavelength 
method  for  Ky-S  band  case  (15.7-3  GHz).  Upper  and  lower 
error  curves  correspond  to  extreme  errors  assuming  ^ 2 dB 
uncertainty  in  Z and  ^25%  uncertainty  in  k.  The  zero  error 
case  corresponds  to  Ng  « .08  cm"^. 
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more  than  60%  for  the  lower  error  case  for  any  fixed  Z.  Hence  it  would 
appear  from  the  above  results  that  the  modified  dual  wavelength  method 
for  the  K^-S  band  system  configuration  is  feasible  and  the  K^-S  band  case 
should  be  avoided, 

5.3.3  Improvement  of  Measurement  Accuracy  Using  the  Modified 
Dual  Wavelength  Method 

Since  the  35  GHz  frequency  is  most  sensitive  to  uncertainties 
of  the  RHS  of  (5.19),  we  will  utilize  this  frequency  for  the  following 
example.  Re-ejqpress  (5.19)  as  follows, 

1 2 (r)  ) 

(U1S)„  . 10  jrf^j  • 2 y Jr 


where  again  the  subscript  o assumes  the  zero  error  case.  The  (LHS)  of 
(5.33)  corresponds  to  the  zero  error  measured  LHS  of  (5.19).  Also  is 
given  by, 

b 

k = a Z ° 
o 0 0 

For  the  case  in  which  measurement  errors  exist,  we  may  express  (5.33)  as 

(“S),  ■ I”  '^ho  jl^j  * 

where  the  subscript  e denotes  the  error  case.  We  may  relate  the  error 

values  to  the  true  values  Z and  k by  the  relations. 

o o 

Ze  - Z^(l  . 6^) 

^ = •'o^l  ^ «k) 


(5.33) 


(5.34) 


(5.35) 


(5.36) 


(5.37) 


We  note  that  if  the  regression  relation  (5.34)  is  obtained  through  inde- 
pendent measurements  of  k and  Z,  only  the  error  in  k (i.e.,  6j^) , influences 
the  regression  relation.  Substituting  the  above  relations  in  (5.35), 
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In  Table  5.10  are  listed  the  conqputed  errors  (5.42)  for  the 
various  cases  of  the  nominal  ranges  of  10.35,  10.95,  and  11.55  km  in 
Fig.  40.  These  latter  ranges  were  selected  as  the  Z levels  are  higher 
resulting  in  the  difference  term  [(LHS)  - (LHS)  ] dominating.  We  note 
that  the  originally  assumed  error  in  is  approximately  predicted  and 

may  thus  be  removed. 

' Of  course,  the  above  cinalysis  assumes  a zero  error  in  the 

measurement  of  the  LHS  of  (5.19).  As  mentioned  previously  since  this  is 
j a relative  quantity,  it  should  be  measurable  to  within  * 0.5  dB  and  only 

i a minor  perturbation  of  the  calculated  value  of  6,  should  result. 

t 

With  the  error  due  to  k reduced,  the  drop  size  distributions 
obtained  via  the  dual  wavelength  method  may  be  recomputed  and  used  to 
establish  regression  relationships  k = aZ*’  at  other  wavelengths. 
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Table  5.10 
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